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FOREWORD 

This  report  documents  the  first  phase  of  the  effort  performed  for  the  Naval  Sea 
Systems  Command  Code  06H1  on  contract  N00024-77-C-6251.  under  Project  Ser- 
ial Number  SF  11  121  160,  Task  20929,  anil  covers  the  period  from  October  1977 
to  October  1978.  It  includes  all  of  the  technical  material  found  in  the  three 
— Quarterly  Reports  plus  work  performed  in  the  fourth  quarter,  connective  material, 
the  conclusions  of  first  phase,  and  the  approach  for  the  second  phase  of  the  study. 
This  study  has  as  its  goal  to  examine  the  application  of  adaptive  filters  to  perform 
passive  sonar  bearing  tracking.  The  first  major  steps  toward  this  goal  have  been 
realized  via  a combination  of  statistical  modeling  and  analysis  that  has  been 
verified  via  computer  simulation. 
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I.  INTRODUCTION 

This  is  the  final  report  on  the  first  phase  of  a study  contract  to  examine  the  applicability 
of  adaptive  filtering  to  passive  sonar  bearing  tracking.  The  study  was  motivated  by 
certain  limitations  in  existing  split-beam  trackers  that  might  be  removed  by  using  an 
adaptive  filter  canceller.  In  particular  1)  an  adaptive  filter  does  not  require  a-priori 
power  spectral  information  on  the  signal  and  noise  fields,  2)  because  the  filter  charac- 
teristics are  adjusted  iteratively  (time-varying),  the  filter  has  the  capability  to  track 
dynamic  (non-stationary)  inputs  and  3)  since  all  the  correlation  information  between 
split  array  outputs  is  contained  in  the  adaptive  filter  canceller  weights,  the  potential 
exists  to  perform  both  broadband  and  narrowband  tracking  simultaneously,  using  all 
the  signal  energy  rather  than  separating  the  energy  as  in  current  tracking  systems. 

This  report  contains  the  results  reported  in  the  first  three  quarterly  reports  with  some 
modifications  and  new  results  obtained  during  the  fourth  quarter  effort.  The  body  of 
the  report  is  organized  so  as  to  present  a summary  of  the  results  of  the  analytical 
effort.  The  details  of  the  analyses  are  found  in  the  Appendices. 

H.  STATEMENT  OF  PROBLEM 

The  work  effort  has  been  focussed  on  studying  the  following  set  of  analysis  tasks  asso- 
ciated with  using  the  adaptive  filter  as  a bearing  tracker: 

1.  0 Analysis  of  the  application  of  adaptive  fillers  to  bearing  tracking. 

1.1  Determine  the  best  method  for  obtaining  the  delay  parameter  t from  the 
adaptive  filter  weights.  Compare  the  bias  and  variance  from  this  method  of 
estimating  r to  that  obtained  from  existing  split  beam  trackers. 

1.2  Determine  the  sensitivity  of  tracking  performance  to  unknown  signal  and 
noise  power  spectra. 

1.3  Evaluate  the  bias  caused  by  plane  wave  interference  on  tracker  performance. 

1.4  Investigate  the  design  of  an  adaptive  tracker  which  can  use  both  broadband 
and  narrowband  sigr  a components  simultaneously  to  estimate  target  track. 
Determine  if  the  ambiguities  in  the  adaptive  weights  can  be  resolved  for 
narrowband  inputs. 

2.0  Analysis  of  adaptive  filters  with  non-stationary  inputs. 

2.1  Determine  the  closed  form  transient  response  of  an  adaptive  line  enhancer 
having  linearly  changing  input  signal  frequency  from  the  matrix  eigenvalues. 
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2.2  Model  broadband  dynamic  Inputs  into  the  adaptive  tracker  and  analyze 
the  time  response  and  mean  squared-error  of  the  cancelled  (error) 
output. 

3.0  Simulation,  test  and  evaluation  using  computer  generated  data. 

•3.1  Model  the  adaptive  tracker  on  the  RAC  simulation  facilities  and  determine 

performance  using  synthetic  data.  Determine  tracker  bias  and  variance  as  a 
function  of  signal  bearing  rate,  S/N  ratio,  and  plane  wave  interference. 
Broadband  and  narrowband  signals  as  well  as  combined  broadband  and  narrow- 
band  signals  are  to  be  used. 


ID.  SUMMARY  OF  RESULTS 

The  results  of  this  study  are  divided  into  seven  parts.  The  first  part  deals  with  the 
selection  of  the  adaptive  filter  configuration  to  perform  bearing  estimation.  Both  sum 
beam  and  split-beam  tracker  configurations  were  studied. 

The  second  part  deals  with  the  extraction  of  the  bearing  (or  delay)  estimate  from  the 
adaptive  filter,  keeping  in  mind  that  the  AF  is  basically  operating  as  a waveform  esti- 
mator. For  broadband  inputs,  the  delay  estimate  is  chosen  to  correspond  to  the  peak 
value  of  the  weights.  The  third  part  deals  with  a frequency  domain  equivalent  model 
of  the  time  domain  adaptive  filter  that  enables  the  variance  of  the  adaptive  filter  delay 
estimator  to  be  calculated.  The  variance  of  the  delay  estimator  is  compared  to  the 
Cramer-Rao  Lower  Bound  (CRLB)  for  the  variance  of  any  unbiased  estimator.  The 
fourth  part  deals  with  practical  implementation  of  the  delay  estimator  that,  because  of 
the  discrete-time  locations  of  the  filter  taps,  differs  from  the  theoretical  continuous- 
time model  whose  performance  was  analyzed  in  Part  3.  The  discrete  time  implemen- 
tation consists  of  an  Interpolation  routine  following  the  adaptive  filter  that  uses  the 
weight  values  to  Interpolate  delays  between  the  taps.  The  degradation  of  the  interpola- 
tion model,  as  compared  to  the  continuous  estimator,  is  investigated. 

The  fifth  part  deals  with  the  performance  of  the  adaptive  filter  bearing  estimator  as  a 
tracker  for  signals  that  are  changing  bearing.  Three  signal  models  are  considered 
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single  frequency  signals  (narrowband  inputs),  broadband  spectrally  white  signals,  and 
broadband  non-white  signals.  In  each  case,  it  is  shown  that  the  adaptive  filter  is  capa- 
ble of  tracking  changing  bearing.  A related  dynamical  problem  is  also  investigated; 
the- response  of  the  adaptive  filter  to  a linearly-varying-frequency  sine  wave. 

The  sixth  part  deals  with  the  simulation  of  the  adaptive  filter  bearing  estimator  in  both 
the  time  domain  filter  configuration  and  its  frequency  domain  equivalent.  For  a wide 
variety  of  simulation  parameters,  the  two  configurations  are  shown  to  perform  equally 
well  and  to  yield  close  agreement  with  the  theoretical  predictions.  It  also  compares 
the  adaptive  tracker  performance  to  that  of  a conventional  tracking  system. 

The  seventh  part  deals  with  the  conclusions  of  the  study  that  are  based  upon  the  results 
of  the  previous  six  parts.  These  conclusions  are  related  to  the  future  tasks  to  be  per- 
formed in  Phase  2 of  the  overall  study. 

Part  I - Selection  of  the  Adaptive  Filter  Configuration  for  a Bearing  Estimator 

It  would  have  been  highly  desirable  to  apply  some  type  of  optimization  theory  to  the 
selection  of  the  structure  containing  the  adaptive  filter  and  constituting  the  bearing 
estimator.  For  example,  what  structure,  using  the  adaptive  filter,  yields  the  minimum 
variance  unbiased  estimate  of  bearing  in  the  class  of  allowable  structures.  Unfortu- 
nately, no  such  optimization  theory  exists  nor  does  it  appear  feasible  to  attempt  to 
derive  a theory.  Instead,  based  upon  some  reasonable  assumptions  about  what  struc- 
tures might  be  useful  for  bearing  estimation,  two  bearing  estimator  structures  were 
selected  on  an  ad  hoc  basis  and  studied  in  detail. 

The  two  bearing  estimator  structures,  sum  beam  and  split-beam,  are  shown  in  Fig- 
ures 1 and  2.  It  is  shown  in  Appendix  I that  the  sum -beam  system  tends  to  obscure 
the  signal  correlation  properties  necessary  for  the  adaptive  filter  to  track.  The  hydro- 
phone noises  produce  residual  correlations  on  the  beams  which  the  adaptive  filter 
treats  as  a signal  process.  To  decorrelate  the  noises,  one  must  form  beams  from 
separate  elements,  so  that  the  two  inputs  to  the  adaptive  canceller  are  not  linear  com- 
binations of  the  same  set  of  processes.  This  motivated  the  consideration  of  the  split- 
beam  processor  of  Figure  2.  It  is  shown  in  Appendix  I that  the  split-beam  system 
offers  the  capability  of  correctly  estimating  the  delay  between  the  two  split-beam  inputs. 
Furthermore,  in  Appendix  IV,  it  is  shown  that  the  performance  of  the  split-beam 
adaptive  bearing  tracker  is  very  near  to  the  theoretical  lower  bound  on  performance  of 
any  tracker.  This  latter  result  lends  strong  support  to  the  ad  hoc  selection  process. 
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Part  II  - Estimation  of  Bearing  Estimates  from  the  Adaptive  Filter  Weights 

The  structure  for  estimating  bearing  (delay)  is  shown  in  Figure  3.  The  left  half-beam 
output  is  selected  to  be  the  desired  signal  and  the  right  half-beam  is  selected  to  be  the 
input  to  the  adaptive  filter.  Ideally  the  filter  should  insert  a delay  equal  to  the  signal 
arrival  time  difference  between  the  two  half  array  phase  centers.  Hence,  one  weight 
(at  the  correct  delay  on  the  delay  line)  will  be  unity  and  all  other  weights  will  be  zero. 
However,  because  the  filter  is  a tapped  delay  line,  only  discrete  values  of  the  time 
delay  can  be  estimated  exactly  by  observation  of  a single  non- zero  weight.  Thus,  de- 
pending  on  the  spectra  of  XL(t)  and  X^(t),  more  than  one  tap  weight  will  be  non- zero. 

If  we  also  include  the  effects  of  noise,  it  is  conceivable  that  several  weights  will  be 
non- zero  even  if  the  true  delay  occurs  precisely  at  a tap  on  the  delay  line.  Hence, 
some  kind  of  algorithm  is  needed  for  determining  the  correct  location  of  the  peak  of 
the  weights  from  noisy  sample  values  of  the  weights  taken  at  discrete  times.  The 
particular  algorithm  is  discussed  in  Appendix  V. 

A theoretical  analysis  of  the  performance  of  the  bearing  estimation  scheme  in  Figure  3 
requires  statistical  knowledge  of  the  filter  weights  after  convergence.  Although  the 
mean  values  of  the  weights  can  be  determined  as  a function  of  the  signal  and  noise 
powers  and  the  algorithm  gain  coefficient  p,  the  variance  of  the  time  domain  weights 
in  the  presence  of  correlated  inputs  to  the  two  half-arrays  cannot  be  determined 
'directly  at  this  time.  Knowledge  of  the  weight  variance  is  necessary  in  order  to  com- 
pute the  variance  of  the  bearing  estimate.  Because  of  this  lack  of  statistical  knowledge 
about  the  system  Figure  3,  an  equivalent  frequency-domain  adaptive  bearing  estimator 
is  studied  in  the  next  section.  The  alternate  configuration  is  analyzed  and  used  to 
predict  the  performance  of  the  system  in  Figure  3.  Subsequent  analysis  and  simula- 
tions verify  the  equivalence  of  the  two  systems. 

Part  III  - Frequency  Domain  Equivalent  of  the  Time  Domain  Adaptive  Filter  Bearing 
Estimator 

Figure  4 shows  an  alternate  structure  for  the  split-beam  adaptive  bearing  estimator 
that  displays  more  clearly  the  ability  of  the  device  to  track  both  broadband  and  narrow- 
band signals.  Most  conventional  trackers  are  designed  to  handle  either  narrowband 
signals  or  broadband  signals  but  not  both.  For  example,  a bearing  tracker  that  uses 
a phase-locked  loop  (PLL)  to  measure  the  difference  in  phase  between  two  half-arrays 
can  effectively  track  narrowband  signals  where  the  phase  of  a carrier  with  slowly 
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varying  envelope  can  be  followed.  However,  broadband  signals,  which  do  not  have 
a well-defined  carrier  are  not  tracked  well  by  a PLL. 


The  split-beam  adaptive  bearing  tracker  can  track  either  broadband  or  narrowband 
signals  since  the  adaptive  filter  does  not  need  to  know  the  signal  statistics  a priori. 

The  delay  of  a broadband  signal  will  be  seen  as  a relatively  large  peak  in  the  weight 
values  at  the  correct  delay  setting.  The  delay  of  a narrowband  signal  will  be  seen  as 
a phase  delayed  set  of  sinusoidal  filter  weights.  In  thikqase,  the  delay  will  display 
itself  most  easily  by  FFT'ing  the  weights  and  selecting  the  largest  output  of  the  phase 
estimator  in  Figure  4. 

The  significance  of  the  above  arguments  is  that  the  split-beam  adaptive  bearing  tracker 
automatically  configures  itself  as  a broadband  or  narrowband  tracker  as  the  situation 
requires.  This  means  that  the  device  is  able  to  track  on  all  the  incoming  signal  energy, 
whether  broadband  or  narrowband.  The  processor  of  Figure  4 is  a method  for  exploit- 
ing this  behavior.  The  peak  in  the  adaptive  filter  weights,  or  in  the  Fourier  transform 
of  the  weights  is  used  to  determine  the  inter-array  delay  or  phase  shift. 

Figure  4 indicates  the  structure  for  extracting  bearing  estimates  from  the  time-domain 

adaptive  filter.  For  comparison,  Figure  5 shows  a frequency  domain  analog  of  the 

time  domain  adaptive  filter.  In  Figure  4,  it  is  noted  that  the  entire  broadband  input  is 

the  input  to  M tap  adaptive  filter.  The  error  signal  over  the  total  band  controls  the 

weight  adjustment  algorithm.  On  the  other  hand  in  Figure  5,  the  input  is  divided  into 

narrow  frequency  bins  with  each  narrowband  input  and  error  signal  controlling  the 

adjustment  of  a single  complex  weight  for  each  bin.  Since  the  weight  adjustment  in 

each  frequency  bin  occurs  independently  of  errors  in  other  frequency  bins,  it  is  not 

obvious  that  the  two  implementations  are  equivalent.  In  Appendix  n,  for  the  case  of 

a broadband  signal  in  noise,  it  is  shown  that  the  implementations  in  Figures  4 and  5 

converge  to  essentially  the  same  steady-state  weight  vector  and  mean-square  error 

if  the  FFT  time  window  in  Figure  5 is  large  in  comparison  to  the  time  delay  difference 

between  the  signals  at  the  adaptive  tracker  inputs.  Furthermore,  by  examining  the 

difference  equations  for  the  mean  value  of  the  weights,  it  is  shown  in  Appendix  II  that 

the  time  responses  are  equal.  The  weight  fluctuations  are  also  compared  and  shown 

2 

to  differ  by  a ratio  that  is  approximately  unity  for  n,Tn  small. 
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Having  verified  the  equivalence  of  the  time  and  frequency  domain  adaptive  bearing 
estimator  configurations,  we  now  proceed  to  determine  the  performance  of  the  fre- 
quency domain  bearing  estimator. 

The  performance  of  the  broadband  adaptive  filter  tracker  requires  the  mean  and  vari- 
ance of  the  weights  with  both  signal  and  noise  present.  The  configuration  initially 
proposed  operates  entirely  In  the  time  domain  using  a tapped  delay  line  (Figure  4). 

This  is  contrasted  with  the  structure  in  Figure  5 wherein  the  beam  outputs  are  pre- 
filtered  using,  for  example  a Fast  Fourier  Transform  (FFT)  and  a single  complex 
adaptive  filter  weighting  is  used  on  each  FFT  bin  on  each  split  beam.  The  time  delay 
between  split  array  phase  centers  corresponds  to  a different  phase  shift  at  each  fre- 
quency. The  adaptively  filter  outputs  can  subsequently  be  recombined  to  produce  a 
broadband  output  if  desired.  The  key  assumption  being  employed  is  that  the  obser- 
vation time  is  relatively  long  compared  to  the  inverse  of  the  bandwidths  so  that  Fourier 
coefficients  provide  an  adequate  second-order  statistical  realization  of  the  broadband 
process  itself. 


If  the  input  processes  are  wide-sense  stationary  over  the  observation  time,  then 
disjoint  spectral  outputs  are  uncorrelated.  Since  the  input  data  is  assumed  to  be  a 
narrowband  gaussian  random  process,  the  disjoint  bins  therefore  provide  statistically 
independent  outputs.  Thus  each  complex  tap  is  operating  on  independent  data  just  as 
the  time  domain  taps  along  the  delay  line  are  assumed  to  operate  on  samples  that  are 
independent  in  time.  This  interpretation  casts  the  analysis  in  a form  that  is  tractable. 
In  Appendix  HI,  the  mean  and  variance  of  the  complex  weight  is  calculated  with  both 
signal  and  noise  present.  This  result  is  then  used  to  obtain  steady-state  tracker 
performance. 


The  principle  results  of  Appendix  m are 
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a - normalized  signal  correlation  between  two  half  arrays 
o 

- signal  power 
u 2 - noise  power 

H algorithm  feedback  coefficient 
E = denotes  statistical  expectation 

Furthermore,  the  steady-state  mean-square  error  is  given  by 


The  above  results  can  be  used  to  determine  the  performance  of  a continuous  version  of 
the  adaptive  bearing  estimator  shown  in  Figure  3.  Figure  6 differs  from  Figure  3 in 
that  the  tapped  delay  line  in  the  adaptive  filter  is  assumed  to  be  infinitely  dense  (i.  e. , 
the  tap  separations  are  infinitesimal).  This  allows  us  to  view  the  weight  values  as  a 
function  of  a continuous  parameter,  say  r,  rather  than  as  a function  of  a discrete 
parameter,  m,  as  in  Figure  3.  The  advantage  of  this  viewpoint  will  be  evident  shortly. 

We  furthermore  assume  the  time  domain-frequency  domain  equivalence  of  Fig- 
ures 4 and  5 is  valid  as  verified  in  Appendix  n and  apply  the  results  of  Appendix  in  to 
Figure  6. 

The  model  and  statistical  weight  analysis  of  Appendices  II  and  III  are  used  to  make  a 
first-order  approximation  analysis  of  the  time-domain  adaptive  filter  when  configured 
as  a bearing  estimator.  This  analysis  uses  the  model  in  Figure  5 for  broadband  inputs. 
The  objective  is  to  operate  the  broad  adaptive  filler  (in  any  implementation,  time  or 
frequency  domain)  and  compare  its  performance  with  the  Cramer-Rao  Lower  Bound 
(CRLB).  The  parametric  behavior  at  low  signal-to-noise  ratios  is  a fundamental 
concern  of  the  analysis.  It  is  shown  in  Appendix  IV  th3t  the  variance  of  the  split-beam 


adaptive  bearing  estimator  is  very  near  to  CRLB  lower  bound  on  the  variance  of  any 
split  beam  estimator. 

The  variance  of  the  estimator  is  within  0.5  db  of  the  CRLB.  Hence,  the. split-beam 
adaptive  bearing  estimator  performs  close  to  the  theoretical  lower  bound  on  the  vari- 
ance of  any  unbiased  estimator.  Hence,  the  adaptive  bearing  estimator  is  nearly 
optimum  and,  in  effect,  makes  a maximum  likelihood  estimate  of  the  bearing  without 
a priori  knowledge  of  the  signal  and  noise  spectra.  The  price  of  unknown  input  spectra 
is  the  0.  5 db  difference  in  performance  in  comparison  to  the  CRLB. 

Part  IV  - Delay  Estimation  with  a Discrete  Adaptive  Tracker 

In  Part  3,  an  expression  for  the  variance  of  the  bearing  estimate  of  a continuous 
adaptive  filter  configured  as  an  adaptive  tracker  was  developed.  The  estimate  is 
based  upon  determination  of  the  peak  of  the  continuous  impulse  response  of  the 
adaptive  filter.  In  practice,  however,  the  adaptive  filter  is  discrete  in  time  and  of 
finite  length,  and  the  peak  of  the  impulse  response  must  be  determined  by  interpo- 
lating between  the  discrete  sample  points.  The  interpolation  process  increases  the 
variance  of  the  bearing  estimate  with  respect  to  that  of  the  continuous  case. 

AppendLx  V describes  a numerical  method  of  determining  the  peak  of  an  interpo- 
lated impulse  response,  and  develops  an  expression  for  the  variance  of  the  esti- 
mate using  this  method. 
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Part  V - Adaptive  Bearing  Tracker  Performance  for  Dynamic  Signals 

The  previous  analysis  has  considered  the  behavior  of  adaptive  bearing  estimator  for 
signals  at  a fixed  bearing  angle;  that  is,  how  does  the  performance  of  adaptive  bearing 
estimator  compare  to  the  CRLB  for  static  signals.  Although  the  static  analysis  gives 
a good  indication  oi  how  the  estimator  performs,  the  real  case  of  interest- is  in  the 
tracking  mode  where  the  bearing  angle  of  the  target  is  changing.  In  Appendix  VI  the 
time-varying  mean  weights  are  derived  for  a signal  that  is  moving  such  that  the  delay 
between  split-array  phase  centers  is  linearly  changing  with  time.  Three  signal  models 
are  considered  - a narrowband  signal,  a broadband  spectrally  white  signal  and  a 
broadband  non-white  signal.  For  the  single  frequency  case,  the  mean  filter  weights 
correct  for  a frequency  shift  and  phase  precession  due  to  the  angular  movement  of  the 
target  in  bearing.  For  a spectrally  white  process,  the  mean  filter  weights  are  a 
travelling  wave  with  a decaying  exponential  envelope.  Tracking  the  bearing  involves 
estimating  the  delay  location  of  the  leading  edge  of  the  weights.  For  the  broadband  non- 
white signal,  the  mean  weights  behave  in  a similar  manner  to  the  previous  case  except, 
because  of  correlation  between  the  taps,  the  peak  of  the  travelling  wave  is  much  larger 
and  the  decaying  exponential  envelope  is  dependent  on  the  signal  dynamics  and  correla- 
tion. Thus,  tracking  non-white  signals  is  easier  than  tracking  white  signals. 

A subject  related  to  the  tracking  behavior  of  the  adaptive  filter  is  treated  in  Appendix  VII. 
The  transient  behavior  of  the  LMS  adaptive  filter  is  studied  when  the  filter  is  configured 
as  a canceller  operating  in  the  presence  of  a fixed  or  variable  complex  frequency  sine- 
wave  signal  buried  in  white  noise.  For  a fixed  frequency  signal,  the  mean  weights  are 
shown  to  respond  to  signal  more  rapidly  than  to  noise  alone.  For  a chirped  signal,  a 
fixed  parameter  matrix  first  order  difference  equation  is  derived  for  the  mean  weights 
and  a closed  form  steady-state  solution  obtained.  The  transient  response  is  obtained 
as  a function  of  the  eigenvectors  and  eigenvalues  of  the  input  covariance  matrix.  Suf- 
ficient conditions  for  the  stability  of  the  transient  response  are  derived  and  an  upper 
bound  on  the  eigenvalues  obtained.  Finally,  the  mean-square  error  is  evaluated  when 
responding  to  a chirped  signal.  The  gain  coefficient  of  the  LMS  algorithm  is  determined 
that  minimizes  the  mean-square  error  for  chirped  signals  as  a function  of  chirp  rate 
and  signal  and  noise  powers. 

Part  VI  - Simulations 

The  last  part  of  this  report  deals  with  simulations  of  adaptive  bearing  trackers  in  both 
the  time  and  frequency  domain  Implementations  followed  by  time  domain  interpolation. 
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The  simulations  serve  three  purposes: 


1)  Verification  of  the  equivalence  of  the  time  and  frequency  domain  implementa- 
tions of  the  tracker, 

. 2)  Support  for  the  theoretical  model  of  the  frequency  domain  implementation, 

. 3)  Provide  initial  sets  of  data  for  design  of  the  adaptive  bearing  tracker. 

Figures  17  - 46  show  the  delay  estimates  as  a function  of  the  various  system  parame- 
ters indicated  on  the  figures. 

In  order  to  verify  the  results  of  the  previous  sections  and  the  appendices,  the  adaptive 
tracker  has  been  simulated  on  the  computer.  The  simulation  is  configured  to  imple- 
ment both  the  time  and  frequency  domain  trackers,  and  provides  for  broadband  or 
narrowband  signal  inputs  with  a variable  delay  between  the  simulated  array  inputs. 

The  capability  is  included  to  band  limit  the  input  via  a finite  impulse  response  (FIR) 
filter  in  each  input  path.  The  delay  may  also  be  varied  linearly  with  time.  A sin  (x)/x 
interpolator,  as  described  in  Appendix  V,  is  included  to  estimate  delays  when  the 
true  delay  lies  between  discrete  filter  taps.  A phase  estimator,  corresponding  to  the 
narrowband  case,  is  also  included.  Figures  13  and  14  illustrate  the  simulation  for 
time  and  frequency  domain  trackers,  respectively. 

A primary  result  of  the  second  quarterly  report  was  that  the  time  delay  can  be  ex- 
tracted from  the  weight  vector  even  at  low  signal- to-noise  ratios  where  the  reduction 
in  mean  square  error  is  small.  For  a broadband  signal  with  a delay  equal  to  an  inte- 
gral number  of  sampling  intervals,  this  would  be  evidenced  by  the  filter  weight  corre- 
sponding to  the  correct  delay  appearing  larger  than  the  other  filter  taps.  The  ability  to 
distinguish  that  the  correct  tap  is  largest  at  input  signal  to  noise  ratios  as  low  as  -20  dB 
has  been  verified  by  simulation.  Figures  15,  16,  and  17  show'  the  time  domain  weight 
vector  as  a function  of  time  for  one  sample  interval  (0.4167  msec)  delay  between  tracker 
inputs  for  a broadband  signal  with  0,  -10,  and  -20  dB  signal- to-noise  ratios,  respec- 
tively. The  2nd  weight,  corresponding  to  one  unit  of  delay,  clearly  has  the  largest 
mean  value  in  all  cases.  Figures  18,  19,  and  20  show  similar  results  with  the  sig- 
nal delayed  3-sample  intervals  between  inputs. 

The  results  of  the  simulations  using  the  frequency  domain  implementation  verify  the 
analytical  prediction  of  the  steady  state  weight  vector  derived  earlier  in  die  report. 
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Table  1 below  gives  the  value  ox  the  largest  steady  state  time  domain,  mean  weight 
for  broadband  examples  simulated  with  the  frequency  domain  tracker.  The  remaining 
weights  should  have  mean  zero  in  steady  state.  These  results  can  be  predicted  using 
equation  (8)  of  Appendix  n.  The  theoretical  predictions  can  be  compared  to  the  simu- 
lation results  given  in  Tables  3 through  8.  These  show  the  time  domain  weight  vector 
as  a function  of  time  for  each  of  the  six  cases  of  Table  1.  The  value  of  the  weights 
shows  good  agreement  with  the  above  analytical  results  once  the  weights  approach  steady 
state,  as  they  have  in  the  low-er  part  of  each  table.  It  should  be  noted  that  each  2400 
iterations  constitutes  one  second  in  these  simulations.  The  tables  also  indicate  that  it 
is  possible  to  determine  the  correct  delay  from  the  weight  vector  for  these  signal-to- 
noise  ratios,  since  the  correct  weight  is  largest  in  all  cases. 


TABLE  1.  PREDICTED  VALUE  OF  LARGEST  STEADY  STATE  MEAN 
WEIGHT  IN  THE  TIME  DOMAIN 


SNR 

V- 

Delav 

Simulation  Results 
in  Table  No. 

1 

Largest  Steady  State  ' 
Mean  Weight 

0 dB 

2-10 

: 

i 

3 

W2  = 0.  4375 

0 dB 

2'10 

3 

4 

W . = 0.  3125 

4 

; -10  dB 

2'10 

2 

1 

5 

Wg  = 0.  0795 

1 -10  dB 

2'10 

3 

6 

W,  = 0.05G6 

4 

-20  dB 

2-14 

1 

7 

i 

W2  = 0.00867 

-20  dB 

2-14 

3 

8 

W4  = 0.00619 

Using  the  delay  interpolator  described  in  Appendix  V,  a number  of  simulations  with 
broadband  signal  have  been  run.  The  inputs  to  the  adaptive  tracker  are  low  pass  fil- 
tered with  a cutoff  of  800  Hz  and  sampled  at  2400  Hz.  A 16-tap  time  domain  filter  is 
used  to  provide  the  weights  to  the  interpolator,  which  has  its  bandwidth  parameter,  B, 
matched  to  the  filter  input  bandwidth,  800  Hz.  Figures  21  through  23  show  the  delay 
estimate,  t versus  time  for  m = 2-10  and  an  actual  delay  halfway  between  the  Sth  and 
9th  taps  (3.541667  msec),  with  signal-to-noise  ratios  (SNRs)  of  10,  0,  and  ■ i.0  dB. 

The  one  standard  deviation  (It)  limits  for  the  estimate  as  predicted  using  the  results 

of  Appendix  V is  also  shown.  Figures  24  through  26  show  the  same  cases  with  a 

-14 

feedback  coefficient  of  \x-  2 , which  reduces  the  variance  of  the  estimate.  Finally, 


m 


■ 


Figure  27  shows  a -20  dB  SNR  case  with  n = 2 These  cases  are  repeated  for  the 
actual  delay  coincident  with  the  8th  tap,  or  3.3333  msec,  in  Figures  28  through  34. 

The  following  observations  may  be  made  concerning  these  simulation  results; 

(1)  For  static  delays,  the  adaptive  tracker  with  interpolation  between  weights  can 
provide  high  accuracy  delay,  and  hence,  bearing  estimation  for  broadband 
targets.  Table  2 below,  converts  the  delay  estimate  accuracies  obtained  in 
the  simulations  shown  in  Figures  21  through  27  to  broadside  bearing  accu- 
racies for  split  arrays  with  a 7.5  and  75  foot  spacing  between  phase  centers. 


TABLE  2.  BEARING  ACCURACIES  FOR  TYPICAL  ADAPTIVE  TRACKER 
SIMULATIONS  WITH  BROADBAND  INPUTS 


SNR 

M 

Figure  No. 

r.m.s.  Bearing  Error 

7.5  Ft  Baseline 

75  Ft  Baseline 

10 

2'10 

21 

.187° 

1.87° 

10 

2'14 

24 

o 

f>* 

m 

o 

.57° 

0 

o 

rH 

1 

22 

.67° 

6.7° 

0 

2"14 

25 

.185° 

1.85° 

-10 

2-14 

26 

.877° 

8.77° 

-20 

2-l« 

27 

5.88° 

58.8° 

As  pointed  out  in  AppendLx  V,  p must  be  set  sufficiently  small  to  assure  that 
the  mean  of  the  peak  weight  vector  exceeds  the  standard  deviation  if  accu- 
rate estimation  is  to  result. 

(2)  The  results  of  the  simulations  are  in  excellent  agreement  with  the  predictions 
given  in  Appendix  V,  as  long  as  p is  set  sufficiently  small.  This  is  most 
evident  in  Figure  11.  Further,  the  criterion  for  setting  p given  in  that 
appendix,  is  verified  by  the  simulation  results.  For  example,  from  Appen- 
dLx  V,  p = 2 10  is  marginal  for  an  SNR  of  -10  dB,  and  Figure  23  shows  that 
occasionally  the  weight  maximum  appears  in  the  wrong  bin. 

(3)  Selection  of  n sufficiently  small  to  make  the  mean  of  the  weights  exceed  the 
standard  deviation  requires  small  p,  and  hence,  long  time  constants  for 
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operation  with  low  SNR  targets.  This  raises  some  question  as  to  whether  an 
adaptive  filter  with  these  long  time  constants  can  track  a moving  target.  It 
can  be  seen  from  the  simulations,  however,  that  estimation  within  the  pre- 
dicted accuracy  is  achieved  in  less  than  a single  time  constant.  For  example, 

- .11 

the  time  constant  for  the  simulations  with  p = 2 is  6.28  seconds,  but  the 

estimate  is  within  predicted  steady  state  limits  within  nominally  1 second. 

This  suggests  an  ability  to  track  dynamics  while  utilizing  large  time  constants 
to  give  high  accuracy  estimates  of  static  targets. 

(4)  Wien  n is  marginal  for  a particular  SNR,  so  that  occasionally  the  peak  weight 
appears  in  the  wrong  tap,  the  estimate  falls  within  predicted  steady  state 
error  limits  most  of  the  time.  In  this  case  large  deviations  from  the  average 
could  be  rejected  as  "wild  points,"  and  the  remaining  values  would  be  good 
estimates. 

In  order  to  relate  these  results  to  existing  trackers,  a first  order  performance  com- 
parison between  the  adaptive  tracker  and  a clipped  two  point  correlator  tracking  scheme 
is  provided  in  Appendix  VUI.  The  latter  is  similar  to  many  existing  tracking  systems. 
The  comparison  is  not  valid  for  high  signal  to  noise  ratio,  and  should  be  used  only  to 
compare  the  two  trackers  in  terms  of  order  of  magnitude  of  performance.  This  is 
because  the  assumptions  necessary  to  make  two  trackers  of  markedly  different  struc- 
tures equivalent  in  some  sense  can  only  be  considered  approximate.  Table  9 shows 
the  bearing  accuracy  predicted  for  both  the  adaptive  tracker  and  the  clipped  two  point 
correlator  for  several  typical  cases,  along  with  the  results  of  the  simulations.  The 
results  are  for  the  same  broadband  signal  considered  in  Table  2,  with  the  target 
stationary  at  broadside  of  a split  array  with  7.5  ft  between  phase  centers.  The  per- 
formance is  comparable  with  the  adaptive  filter  having  slightly  smaller  variances  for 
all  cases. 

Appendix  VI  of  this  report  considered  the  ability  of  the  adaptive  tracker  to  track 
linearly  varying  delays.  Such  a delay  is  representative  of  a target  with  constant  bear- 
ing rate  passing  through  broadside  relative  to  the  split  array  centers.  Simulations 
were  run  with  the  same  adaptive  tracker  interpolator  configuration  used  in  the  static 
runs,  but  with  the  delay  varying  linearly  with  time.  Again,  the  input  is  broadband, 
and  both  tracker  inputs  are  bandlimited  to  800  Hz.  The  adaptive  filter  has  16  taps, 
and  the  interpolator  bandwidth  is  matched  to  that  of  the  input. 
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TABLE  9.  COMPARISON  OF  BEARING  ACCURACY  OF  AN  ADAPTIVE 
TRACKER  WITH  CLIPPED  TWO  POINT  CORRELATOR 
TRACKER  WITH  EQUAL  TIME  CONSTANTS 


* 


SNR 

r.m.s.  Bearing  Error 

Predicted, 
adaptive  tracker 

Simulation, 
adaptive  tracker 

Predicted,  clipped 
two  point  correlator 

0 

o-lO 

.819° 



2.57° 

0 

2"14 

.202° 

. 185° 

.642° 

-10 

to 

1 

►— 

1.3° 

.877° 

3.54° 

-20 

2-16 

5.8° 

5.88° 

16.25° 

J 

In  Appendix  VI,  it  was  shown  that  the  mean  weight  vector,  as  a function  of  time  is  a 
peak  at  the  actual  instantaneous  delay  between  array  halves,  moving  through  the  filter 
at  the  delay  rate  of  change,  and  with  an  exponentially  decaying  trailing  edge.  This 
analysis  is  verified  by  Figure  35,  which  shows  the  weight  vector  at  5000,  0000,  and 
7000  iterations  for  a 10  dB  signal  and  a delay  changing  at  1.0  msec  per  second.  Fig- 
ure 36  shows  another  view  of  the  behavior  of  the  weight  vector  with  time.  Here,  the 
delay  is  changing  at  0.2618  msec  per  second,  p = 2“10,  and  the  signal-to-noisc  ratio 
is  10  dB.  The  peak  of  die  weight  vector  can  be  clearly  seen  moving  through  the  filter 
with  time. 

A more  quantitative  assessment  of  the  tracking  behavior  of  the  adaptive  tracker  can 
be  had  from  Figures  37  through  39,  showing  the  delay  estimate  as  a function  of  time 
for  a linearly  increasing  delay  of  261.8  psec/sec.  This  is  representative  of  a l°/sec 
target  at  the  split  array  outputs  of  an  array  with  75  teet  between  phase  centers.  Fig- 
ures 37  and  38  show  the  estimate  for  the  261.8  psec/sec  case  with  p = 2“ 10  and  SNRs 
of  10  dB  and  0 dB  respectively.  The  same  rate  of  change  is  shown  for  p = 2-^  and 
SNR  of  -10  dB  in  Figures  39. 

Part  VII  - Conclusions  and  Relationship  to  Phase  2 

It  can  be  concluded  that  the  application  of  an  adaptive  filter  to  split  beam  bearing 
tracking  has  the  potential  to  offer  performance  Improvements  as  discussed  in  the 
introduction,  in  the  static  case,  bearing  estimates  are  comparable  to  the  Cramcr- 
Rao  Lower  Bound.  In  the  dynamic  case,  time  averages  of  simulated  weights  behave 
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as  the  theoretical  predictions  for  the  mean  of  the  weights.  In  addition,  both  narrowband 
and  broadband  signals  can  be  tracked  with  the  same  processor,  suggesting  that  all  tar- 
get energy  can  be  used  for  tracking. 

The  results  reported  herein  strongly  motivate  continuing  the  effort  in  the  following 
areas  of  investigation: 

(1)  Analyze  the  sensitivity  to  plane  wave  interference,  correlated  noise  (both 
spatially  and  temporally),  multipath,  unknown  statistics,  and  tracking 
dynamics. 

(2)  Through  analysis  and  simulation  determine  the  effects  of  the  following  on 
the  adaptive  tracker  performance: 

(a)  Multipath 

(b)  An  array  of  sensors 

(c)  Multiple  targets  to  be  tracked  at  different  bearing  rates 

(d)  Multiple  interferences  to  be  removed 

(e)  The  number  of  filter  taps  required 

(f)  The  design  of  filter  feedback  and  gain  control  parameters 

(3)  From  the  above  analysis  derive  a system  configuration  and  simulate  on  a 
computer.  Evaluate  the  system  performance  using  computer  generated  data 
and  actual  sea  data  recordings  to  be  provided  by  the  Naval  Sea  Systems 
Command. 

(4)  Compare  systems  performance  from  the  above  evaluation  with  the 
Cramer-Rao  Lower  Bound  and  with  existing  Bearing  Deviation  Indicator 
tracker  performance. 

(5)  Formulate  an  optimal  tracker  structure  based  on  the  above  evaluation  and 
suggest  Implementation  into  real  time  hardware. 


Figure  2.  Configuration  for  Split-Beam  Adaptive  Tracker 


figure  Split  Hcam  Adaptive  Hearing  Kstimator  Slruetiire 


Figure  4.  Split-Beam  Adaptive  Bearing  Estimator  Structure 
(Alternate  Configuration) 
for  Both  Broadband  and  Narrowband  Signals 
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Figure  7.  Adaptive  Filter  Impulse  Response  for  Exact  Steering 

Vectors  and  Spatially  Uncorrelated  Noise,  P = 0 
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Figure  9.  Steady-Stale  Adaptive  Filter  Weights  as  Time  Changes, 
Tracking  the  Changing  Time  Delay 


Figure  10.  Envelope  ot  the  Tracker  Weights  as  Thev  Move  in  Time 
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Figure  14.  Frequency  Domain  Adaptive  Filter  Simulation  Block  Diagram 


Figure  17.  Time  Domain  Weights  versus  Time 
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Figure  22.  Delay  Estimate  versus  Time  (Actual  Delay  = 3.541667  mSEC) 
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Figure  23.  Delay  Estimate  versus  Time  (Actual  Delay  = 3.541667  mSEC) 


Figure  25.  Delay  Estimate  versus  Time  (Actual  Delay  = 3.541667  mSEC) 


Figure  26.  Delay  Estimate  versus  Time  (Actual  Delay  = 3.541667  mSEC) 
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Figure  27.  Delay  Estimate  versus  Time  (Actual  Delay  = 3.541667  mSEC) 
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Figure  31.  Delay  Estimate  versus  Time  (Actual  Delay  * 3.3333  mSEC) 
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Figure  34.  Delay  Estimate  versus  Time  (Actual  Delay  28  3.3333  mSEC) 
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figure  35.  Weight  Vector  as  a Function  of  Time  (n)  for  Linearly  Varying  Delay  and 
L:  naoband  Input.  M * No.  of  Taps  = 16.  SNR  3 10  dB.  m - 2-^,  C = 1.0  msec/sec 
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APPENDIX  I 


SUMMED  BEAM  VS  SPLIT  BEAM  BEARING  ESTIMATION  . 

The  steady  state  results  for  adaptively  cancelling  one  beamformer  output  from  another 
are  examined  for  two  beams  in  a single  summed  array  and  for  identically  steered 
outputs  of  a split-beam  array.  In  both  cases,  the  filter  weights  and  the  resulting  mean 
square  error  are  derived.  A plane  wave  signal,  a plane  wave  interference  and  noise 
are  present  at  each  hydrophone.  Energy  arrival  angles  are  initially  assumed  to  be 
known  with  one  objective  being  to  assess  the  sensitivity  of  the  process  to  errors  in 
steering  for  use  as  a bearing  tracker.  The  power  spectra  for  all  waveforms  is  kept 
completely  general,  but  is  assumed  to  be  known  a priori. 

Summed  Beam  System 

The  situation  being  studied  is  shown  in  Figure  1.  The  signal,  interference  and  noise 
spectra  at  the  kth  array  element  are  denoted  by  s(u-),  i(u.-),  and  n^a.')  respectively. 

The  waveforms  are  zero  mean  and  mutually  uncorrelated.  At  any  frequency,  = 2-f, 
the  vectors  of  these  waveforms  at  all  the  hydrophones  can  be  expressed  as  follows: 


i(«)  = Djfu:)  i(a>) 


s(w)  = Dg(u>)  s(w) 


I(w)  = 


S(w)  = 


/ 


ej*I  \ 

ej2*I 


\ e^l/ 

I eJ*s  \ 
ej24,s 


\ejm4,S 
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N(w)  = 


/ nl<">  \ 

f n2(w)  I 


n (w> 

\ M 


where 

*1  = u c sin 
d 

© = u — sin 
s c s 

<pj  = interference  energy  arrival  angle 

<p  = signal  energy  arrival  angle 
s 

i(w)  = interference  spectrum  at  frequency  u 
s(oJ)  = signal  spectrum  at  frequency  u> 
d = uniform  element  spacing 

c = sound  speed 


The  progressive  delay  vectors,  Dj(w)  and  Dg(w),  represent  the  propagation  of  the 
plane  wave  interference  and  signal  along  the  uniformly  spaced  M element  line  array. 


Two  summed  beams  are  formed  using  uniform  shading.  The  steering  vectors  are 
D^uj)  and  D2(cj),  where 


D,(«) 


7ej*£  \ 

ej2<t>0 


\ejM^  / 


, 2=1,2 


The  beam  steered  in  the  neighborhood  of  the  signal  is  the  target  beam,  denoted  by 
d(w).  The  reference  beam  steered  to  the  interference  is  the  input  to  the  adaptive 
filter,  denoted  by  x(w). 
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d(w)  = D^+  (u>)  Dj(oj)  i(w)  + D*(u)  Dg(u>)  s(u)  + (u)  N(w) 

x(u)  = D2+(w)  Dj(w)  i(u/)  + D2+(u)  Ds(u)  s(w)  + D 2 (w)  N(u) 

where  the  symbol  + denotes  the  complex-conjugate  transpose  operation.  To  form  the 
steady  state  complex  weight  at  a single  frequency  requires  calculating  the  cross- 
spectrum of  d with  x and  the  power  spectrum  of  x.  Thus, 

E [d(w)  x*(w)l  = D + [Dj  Dj  * Pj  + DgDs"  Ps  + QPn]  D, 

E [x(«)  x*(u)]  = D2+[DiDi+Pi  + DsDg'Ps  + QPn]  D2 

where  the  frequency  arguments  on  the  D's  have  been  dropped  fo~  convenience,  and 
where 

E[NN*j  = QPn 

Q = normalized  spatial  covariance  matrix  of  the  noise 
The  complex  weight,  Bfa;),  is 


B(u) 


E[dx*J 

E[xx*] 

Dl^P{Pl  + DsDs\*^  P2 

Vl°]DrpI  + D.Ds*Ps  + «PniD2 


The  error  e,  is 


e (co)  = d(u>)  - B(uj)  x(u>) 

= [ D + - BD*]  [Dj  i(uj)  + Dg  s(w)  + N(w)J 
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and  the  mean  square  error,  E[  |e|  ],  is 


E[  |e  | 2]  = (Dx+-  BD2+)  [DjDj'Pj  + DgDs  Ps  + QPJ  (D^-B~D2)  . 
To  facilitate  writing  E[  |e  j —J , let  us  define  the  following  spatial  responses. 


D1  DI 

V>I 

crDs 

°2\ 

Dl'D2 


= response  of  array  steered  to  v , to  energy  arriving  at  ^ 

= response  of  array  steered  to  <?2,  to  energy  arriving  at  <f> j 

= response  of  array  steered  to  , to  energy  arriving  at  <P 

1 s 

= response  of  array  steered  to  <p0,  to  energy  arriving  at  <P 

S 

= response  of  array  steered  to  <P^,  to  energy  arriving  at  <?2 


For  a uniformly- shaded,  uniformly-spaced  line  array. 


“s  = e 


• -M+l,  . , , , . M . . , . 

(—)  (\~%)  sm  j - 6g) 

sin  1/2  (4^-4^) 


• ,M+1.  , . Sin  2 <*2  ~ »«> 

p = e"J  ( 2 } (<t>2  ~<bs)  sin  1/2  (d>2  -<y 


-j  (®i  -<ty  s*f(W 

sin  1/2  (4>j-4>2) 


For  a j and  (3 . substitute  4>x  for  4>  • 
I I Is 


The  mean  square  error  can  then  be  written  as 


E[  | e | 2]  = loj 1 2 Pz  + lasl2Ps  + D1TQD1Pn  + lB|2|^(PI)2PI  + |Psi2Pg 

T | + * - _ } 

* °2  QD2P»  J’2  Re  | B PI  +“s?s  Ps  + D1  ®D2Pn>  | 


T 


where 

3/i^iVAVWii. 

wWvw, 

For  the  special  case  where  the  noise  is  uncorrelated  spatially  from  element-to-element, 
the  matrix  Q=I,  the  MxM  identity,  and  the  mean-square  error  is: 


! <M2-G122>P„2  + PSPI  <GlsG2I  ' G1IG2s>2 
+ P,P„  <MGls2  + MG2S2  - 2 GlsG2SG12> 


E[|e \A) 


* PlP,  <MQ2I2  * MG1I2  * 2 G1IG21G12> 
G2I2  PI  * G2S2  Ps  + MP„ 


where 


sin  | 


Jab  ” sin  1/2  (*a™<^)) 


and  the  complex  weight,  B(u),  is 


B(u)  = G1»V. * ,-iirMtrti 

°2I  PI  + G2s  Ps  + ^Pn 


An  interesting  example  occurs  for  the  case  of  correct  steering,  for  which  <t>  = 
and  4>2  = . For  this  case 


Gls  " G2!  " M 
G1I  = G2s  = G12 
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and 


E[U|2]  = 


(M2-G122)  ip/  + PSP,  (M2-G122)  t MP,Pn  ♦ MPj  Pn]. 


M^Pt 


G12  Ps  + MPn 


Of  additional  interest  is  the  case  of  a single  plane  wave  broadband  interference  to  be 

cancelled.  This  can  be  evaluated  by  setting  P =0  and  examining  the  extent  to  which  P 

s 

is  removed  from  the  mean  square  error.  Thus 

-j  (^y1)  <VV  g 

B(oj)  = e * 1 ^ _12. 

M 


I 


and 


Note  that  for  this  case  the  interference  is  completely  removed  and  the  weight  is 
independent  of  the  interference  power.  The  reason  for  this  is  that  the  noise  at  the  two 
beam  outputs  has  a deterministic  correlation  placed  on  it  by  the  beamformer.  This 
correlation  is  exactly  the  same  as  the  correlation  on  the  interference,  which  is  the 
sidelobe  gain.  Thus,  the  filter  weights  due  to  the  noise  correlation  are  the  same  as 
the  weights  needed  to  cancel  the  interference.  When  interference  is  added  to  the 
acoustic  field,  it  is  cancelled  regardless  of  the  interference-to-noise  ratio,  because 
the  filter  weights  do  not  have  to  change  to  accommodate  the  interference  correlation. 
This  will  not  be  the  case  if  the  noises  are  correlated  spatially  or  if  the  steering 
vectors  are  not  precisely  aligned  to  the  energy  arrival  angles.  Then  the  noise  corre- 
lation on  the  beams  are  not  exactly  the  same  as  the  interference  correlations  and  the 
adaptive  filter  will  have  to  change  to  cancel  the  interference.  For  small  interference- 
to-noise  ratios,  the  filter  would  not  change  to  cancel  the  interference  because  the 
resulting  increase  in  noise  power  at  the  output  would  exceed  the  reduction  in  output 
power  due  to  removing  the  interference. 
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The  extent  of  the  filter  weights  in  time  is  twice  the  maximum  time  difference  for 
propagation  from  the  reference  beam  to  the  target  beam.  For  example,  when  the  two 
beams  coincide,  i.  e. , sin  = sin  <p2,  there  exists  just  one  weight  at  t = 0.  If  the 
reference  is  at  forward  endfire  and  the  target  is  at  aft  endire,  then  sin  P^-sin^  = 2 
and  the  total  filter  delay  is  2M^  , which  is  twice  the  time  for  the  wavefront  to 
propagate  across  the  aperture  of  the  array. 

It  can  be  concluded  that  the  correlation  on  the  beams  due  to  the  noise  is  helping  to 
cancel  the  interference.  This  is  a desirable  result  when  performing  a post-beamformer 
interference  cancellation  (PIC),  but  it  obscures  the  ability  to  identify  signal  correla- 
tion properties  for  tracking.  The  noises  produce  residual  correlation  on  the  beams, 
which  the  adaptive  filter  treats  as  a signal  process.  To  decorrelate  the  noises  one 
must  form  the  beams  from  separate  elements,  so  that  the  two  inputs  to  the  adaptive 
canceller  are  not  linear  combinations  of  the  same  set  of  processes.  This  motivates 
considering  the  split-beam  processor,  in  the  next  section. 

S^jlit-Beam  System 

The  split-beam  configuration  is  shown  in  Figure  2.  The  2M-array  elements  are 
divided  equally  and  a beam  is  formed  with  both  sub-arrays  steered  to  the  target.  The 
objective  is  to  subsequently  process  to  get  a refined  estimate  of  the  signal  energy 
arrival  angle  and  to  automatically  follow  bearing  and  spectral  changes. 
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n 


The  adaptive  filter  steady  state  impulse  response  is  the  inverse  transform  of  the 
single  frequency  complex  weights.  The  impulse  response,  which  is  approximated  by 
the  discrete  filter  weights,  is 


fm**j*nt 


. M , , d , . 


, ra  ..M+l.  d . sin  -s  w—  (sin% -sin<p0)  . . , 

_ _L  / ) wj  (sin  *x- sin  *2)  2 c v 1 2;  ^jut  dw 

sin  1/2  w—  (sin<a-sin<K>0)  ^ 

-to  C 1 <2 

w M 

’Til  ^ I • ’"VV k r. 


_ 1 V / jwt  -jw£  (sin?  - sin<p„)  _dw 

■ 1 Z / 2TT 


k=l  ' 


M 

= 1 V 


^ 2, 6 c*  ■ *4  <sin’5i  ■ sin<p2)) 

k=l 


The  filter  weights  are  shown  in  Figure  7 for  this  special  case. 
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The  adaptive  filter  will  converge  to  minimize  the  mean  square  error,  i.  e.  to  minimize 
the  power  in  e . This  implies  that  implicit  within  the  filter  is  estimation  of  the  trans- 
formation to  make  the  waveform  x look  as  much  as  possible  like  d.  If  propagation 
delays  or  those  shifts  are  the  only  differences  between  the  half-array  outputs,  then 
one  should  be  able  to  extract  this  information  from  the  filter  weights  and  the  filter 
should  be  able  to  track  reasonable  dynamics. 


Analysis  of  the  steady-state  performance  in  Figure  2 procedes  in  the  same  way  as 
the  previous  section.  The  signal  energy  arrival  angle  is  vji  , from  which  a vector  of 
propagation  delays,  Dg(u;),  is  formed.  The  steering  direction,  from  which  the  vector 
D(gj)  is  formed,  is  close  to  but  not  exactly  Yg.  Thus 


D,(">  - 


jui-  sin  4* 
eJ  c s 

jw  2-  sin>v 
e c s 


JwMp  sink's 


d 

ju;  - sin  9 
e c 


, D(w)  = 


jojM  — sin 


and  just  as  before, 

d(u)  = D (u>)  Dg(w)  s(u>)  + D 

x(w)  = D (u)  e^*s  s(uj)  + D^N  (cj) 

S X 


r 


where  Nj  and  are  uncorrelated  and  all  processes  are  zero  mean.  To  form  the 
complex  filter  weights  one  needs  the  cross  correlation  of  d with  x and  the  auto- 
correlation of  x. 


E (d(u)  x*(w)J  = D+D  D * D e"-*1^  S(u) 
s s 


E [x(w)  x*(w)]  = DTD  D DS(w)  + D QDS  (w) 
s s S II 


where  Sg(w)  and  Sn(w)  are  the  power  spectra  of  s(w)  and  Nx(u;),  and  where  for  con- 
venience the  frequency  argument  of  the  ste ar in g vectors  has  been  dropped.  The 
complex  weight  is 

b,u)  = ay *•  <->3  . -*D^  * jM*Sss(“) 

' ’ E[x(u)  x*  (w)] 


D D. 


Ss(w)  + D~  QD  SQ(w) 


For  the  case  where  the  noise  is  uncorrelated  spatially  from  element-to-element,  the 

matrix  Q=L  Also,  let  a (y  denote  the  spatial  response  of  the  array  steered  to  v to 

energy  arriving  at  u . Thus 
s 


ct  (ui)»  D (oj)  D (u>) 
s 


and 


j 

2 o - j M u>-  sin 


l°l  e 


B(w)  = 


s Sg(«) 


! 2 


°|  Ss  (u>)  + MSn(u) 


The  resulting  error  power  spectrum  is 


E[  I € (UJ)  |2]  = MSn(w) 


2|a  (ui)|  2 Ss(cu)  + MSnM 


o (u)  I Ss(w)  + MSn(w) 


The  filter  transfer  function  , h(t),  is  the  inverse  Fourier  transform  of  B(w).  It  is 
this  function  which  the  adaptive  filter  attempts  to  approximate  with  its  finite  taps. 
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h(t) 


r 


)utD,  . du 
e 1 B(w)  ^ 


-J 


'*  e Jw(t-M  £ sin  vs) 


Sg(W)  |o  (W) 


|«(u)  I ss(w)+MSn(w) 


dw 
2 r 


Li  the  broadband  case  both  S (w)  and  S (u)  can  initially  be  considered  to  be  spectrally 

s ** 

flat,  i.  e.  independent  of  u.  Then 


h(t) 


-f 


ju(t  - sin'+'g)  .la}*  ^ ' 


a(w)|-Ps+MPn 


dcu 

2tt 


Notice,  that  even  with  white  input  spectra,  the  power  spectrum  at  the  beamformer 
output  is  colored  by  the  frequency  dependence  of  the  spatial  response,  a (w).  Since 

_j  (Mli)  (sin+-sin^s)  sin  y Wy(sin9-  sinvg) 

a (a>)  — s 2 wd  | . 

sm  2^-(sin^-  sinog) 

for  a uniformly  spaced  and  uniformly  shaded  line  array,  it  can  be  seen  that  only  for 
perfect  steering,  i.  e.9  = 9g  does  o(uj)  not  depend  on  frequency.  For  that  special 
case. 


M2P„ 


h(t) 


M(MPg+P 


J 


jw(t-  Md  sin  9 J 


da> 

2 - 


o = 9 


M2P 


M^P  +MP 
s n 


rrt  Md  . . , 
6(t  - — sin9g) 


MPg/P, 

1+MP"  /P 
s n 


r 4-“si"0 


and  the  adaptive  filter  forms  a single  weight  at  the  correct  inter-phase  center  delay, 
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APPENDIX  n 

COMPARISON  OF  TIME  AND  FREQUENCY  DOMAIN  ADAPTIVE 
TRACKERS  WITH  BROADBAND  SIGNAL  INPUTS 


1.  Steady  State  Equivalence  of  the  Time  and  Frequency  Domain  Adaptive  Trackers 
with  Broadband  Signal  Inputs 

This  appendix  considers  the  equivalence  of  the  time  domain  and  frequency  domain 
implementations  of  the  adaptive  bearing  estimator,  shown  in  Figures  4 and  5,  respec- 
tively. In  the  case  of  a broadband  signal  in  noise,  it  is  shown  here  that  the  two  imple- 
mentations converge  to  essentially  the  same  steady  state  weight  vector  and  mean 
square  error  if  the  FFT  time  window  is  long  in  comparison  to  the  time  delay  differ- 
ence between  the  signal  at  the  adaptive  tracker  inputs. 

For  the  broadband  signal  case,  the  time  domain  inputs  to  the  tracker  can  be 
written  as 

x(t)  = s(t)  + nx(t) 

d(t)  = s(t  - ATS)  + nd(t)  0) 

where  s(t),  n (t),  and  nd(t)  are  zero  mean,  stationary  random  processes,  uncorrelated 
in  time  and  independent  of  each  other.  Let  the  delay  between  the  two  inputs,  A , be  a 
multiple  of  the  sampling  interval,  Ts.  An  FFT  is  performed  on  both  x(t)  and  d(t)  every 
RTg  seconds,  so  that  each  frequency  domain  filter  iterates  every  RTS  seconds.  With 
M the  FFT  size,  R = M corresponds  to  FFT  processing  without  overlapping  or  gaps 
between  time  windows,  R < M corresponds  to  overlapping  FFTs  and  R > M indicates 
gap  processing. 

The  input  to  the  adaptive  filter  in  the  k—  FFT  bin  is  given  by 


M-l 

z 

nK) 


s [ (n 


iR)  Ts] 


+ n„  [ (n  + iR)  Tj 


. 2r  . 

-J-jT  nk 


n-i 


e 


(2) 


M-l 


D|_(i)  = £ 
n-0 


S t(n  + iR  -A)  Tg]  + nd  [(n  + iR)  Tg]  e 


• 2"  , 
■J  — nk 


(3) 


th 


The  adaptive  filter  in  the  k FFT  bin  has  a single  complex  weight  with  a mean  steady 
state  value  given  by 


Wk'Rdx'/Rxx- 


(4) 


where 


Rxx*  ■ E [V  X'k  «] 
Rdx*  - E [Dk  « \ <*>] 


(5) 


For  the  particular  case  being  considered  here,  the  covariances  are  as  follows. 


xx 


rM-l  M-l 
£ £ 

s[(n  + iR)  Ts]  + n [(n  + iR)  Ts] 

o 

II 

g 

o 

II 

c 

-1 

L J A 

S [(m  + iR)  Tg] 


2 rr 


nx  [(m  + iR)  Ts] 


-]  m 


(n-m)  k 


- M<7^  + Ma  ^ 
s n 


(6) 


using  the  independence  of  s(t|  and  nx(t),  and  the  fact  that  both  are  uncorrelated  in 
time.  Also 


Rdx 


rM-l 

‘=E  j £ 
Ln  = 


M-l  M-l 
*=F  I Z . £ 

o m = o 


nx  [(m  + iR)  Ts] 


S [(n  + iR  -A)  Tg]+  nd  [(n  + iR)  Ts]  S [(m  + iR)  Ts] 
2 it  -|  2 -i  "Tn“^k 

-j~sr(n'm)kJ=(M-A)frs  e 3 


This  results  because  e|  s £(n  + ik)  Ts]  s*  [(m  + iR  +A)  Ts]  | = 0 except  when  m * n -A, 
and  this  occurs  only  (M  -A)  times  over  the  range  of  the  double  summation. 
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Therefore  the  mean  steady  state  weight  in  the  kth  frequency  bin  is 


_ (M  -A)  <7  -j  — Ak 

W = 5 e 

k „ . 2 2 

M (u  + a ) 

' s nx 


The  error  output  for  this  weight  is 


(M  -A)  a 2 -j  M Ak  „ 

‘ k(i)  = Dk (i>  - M -2;-— 2, e \ (i> 

M < *s  + ffnx  > 


so  the  mean  square  error  in  the  k bin  is 


In  order  to  compare  the  steady  state  results  with  that  of  the  time  domain  filter, 

we  need  the  inverse  FFT  of  the  vector  of  complex  frequency  domain  weights.  F rora 

th 

(8),  it  can  be  seen  that  the  n time  domain  weight  will  be 


i — j ~trnk 
j E - 


n * A , 


(M  -A)  a 


n = 0,  1,  ...  M-l 


M (<r  2 + cx  2 ),  n - A, 
s nx 


The  corresponding  total  mean  square  error  is 


' s nx J 


These  are  to  be  compared  to  values  attained  by  the  time  domain  weight  vector, 


0 


n = 0,  1 


n = A , 


M-l 


and  m ean  square  error, 


(13) 


(14) 


Comparing  (11),  (12),  (13),  and  (14),  it  can  be  seen  that  for  A=  0 the  two  filter 
implementations  yield  identical  steady  state  results.  Further,  if  the  delay,  A , is  a 
small  fraction  of  the  FFT  size,  M,  then 


M - A 
M 

and  the  two  implementations  are  essentially  the  same  in  steady  state.  This  places  a 
specification  on  the  size  of  the  FFTs  to  be  used  in  the  frequency  implementation  of  an 
adaptive  tracker.  It  is  also  interesting  to  observe  that  the  steady  state  of  the  frequency 
domain  implementation  of  the  adaptive  tracker  does  not  depend  upon  R.  This  means 
that,  at  least  with  respect  to  steady  state  performance,  use  of  overlapping  or  gap  FFT 
processing  will  not  affect  the  tracker  performance. 
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2.  Transient  Comparison  of  the  Time  and  Frequency  Domain  Implementations 


The  two  implementations  can  be  compared  by  examining  the  difference  equations 
for  the  mean  value  of  the  weights  in  both  cases  using  the  inputs  in  the  previous  section. 
For  the  time  domain,  the  difference  equation  for  the  m^  weight  is 


VV  (n  + I)  = W (n)  + p 
m ' m 


d(n)  -Z  Wp  (n)  Xp  (n) 


X (n) 
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(15) 


which  when  averaged  is 
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From  the  previous  section. 
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where 


( 1 

if  m = A 

= ) 

(19) 

m A 'i 

( 0 
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In  the  frequency  domain,  the  difference  equation  for  the  mean  of  the  complex  weight 
on  the  k1*1  FFT  bin,  is 


E [ WR  (n  + 1)  ] = E 

Wk  (n) 

h-'jt 

= E 

< <n>] 

j 1 - H E 

!♦ 
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These  terms  were  also  evaluated  in  the  previous  section,  where 
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This  expression  can  be  inverse  transformed  to  return  to  a time  domain  weight  vector  for 
comparison  with  the  time  domain  implementation.  By  transforming  back, 


E [ wm  (n.i)]  - E [Wm  <n>]jl  -MM  ( » 4 . »„2>]  - M (M  -4 ) 6m  A (22. 

n-e 


Equations  (19)  and  (22)  are  now  to  be  compared  for  transient  response  and  weight 
fluctuations.  This  is  most  readily  done  for  the  A = 0 case.  If  the  feedback  coefficients 
are  distinguished  from  one  another  by  letting  pT  and  Ppbe  the  time  and  frequency 
domain  coefficients  respectively,  then 


^ p = Mp  (23) 

to  equate  the  two  difference  equations.  The  steady-state  results  than  agree  exactly, 
where  for  A = 0 the  windowing  effects  are  neglected. 

The  time  constant  for  the  adaptive  algorithm  is  T = • The  input  power, 

c ^ in 

Pin,  is  larger  in  the  time  domain  by  a factor  of  M,  which  is  precisely  the  factor  by  which 
Mp  is  larger  than  p^.  Thus  the  p P^n  product  is  the  same  for  both  systems  and  there- 
fore the  time  responses  are  equal. 

The  fluctuations  of  the  weights  can  be  compared  using  results  derived  in  the 
second  quarter  report.  For  the  time  domain  filter  with  uncorrelated  inputs 


lim  Varf  W (n)  I Hoi  = 
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M T 'T“n,  T 
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and  for  the  frequency  domain 
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Again,  for  the  time  domain  implementation,  (Tn2  j is  in  the  full  band,  whereas  for  the 
frequency  domain  vn<“  p is  in  the  FFT  bin,  so  that 


2 2 
<t  * M <x 


Thus  for  the  frequency  domain 
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and  for  the  time  domain 
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The  variances  differ  by  the  ratio 


2-2V" 

2 -(^K2 


which  is  approximately  unity  for  small. 

It  is  thus  concluded  that  for  all  practical  purposes,  the  frequency  domain  adaptive 
filter  model  is  equivalent  to  the  time  domain  adaptive  filter. 


APPENDIX  ID 

MEAN  AND  VARIANCE  OF  A SINGLE  COMPLEX  TAP  WEIGHT 
WITH  SIGNAL  AND  NOISE  INPUTS 

In  narrow  bands  at  the  FFT  output,  the  adaptive  filter  is  a single  complex  weight. 

The  weight  update  equations  are  given  by 

W(n+1)  = W(n)  + K<(n)  x *(n) 

♦ (n)  = d(n)  - y(n)  = error  (cancelled)  waveform 

y(n)  = W(n)  x(n)  = filter  output 

d(n)  = reference  waveform 

x(n)  = input  data  sequence 

W(n)  = complex  weight 

(j.  = feedback  coefficient. 

Thus  the  filter  is  described  here  at  the  nth  time  iteration  on  the  data  sequence  of  a 
particular  FFT  bin.  One  such  device  operates  on  each  of  the  FFT  bins  spanning  the 
signal  bandwidth. 

Consider  again  the  structure  in  Figure  5,  where  the  common  plane  wave  component 
is  denoted  by  I(n).  Then 

d(n)  = al(n)  + n^n) 

x(n)  = I(n)  + n2(n) 

where  a is  the  complex  coefficient  for  the  inter-array  phase  shift. 
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It  is  assumed  that  the  input  sequences  are  zero  mean,  gaussian,  independent  of 
one  another,  and  uncorrelated  in  time. 

The  error  sequence  t (n)  is 

«(n)  = ol(n)  - W(n)I(n)  + n^n)  - W(n)  n2(n) 

The  mean  square  error  is 

E [ | « (n)j 2 J = Ej|[el(n)  - W(n)I(n)]|  2 + [j  *1  (n)j2  + | n2  (n)W  (n)|2  j J 

The  variance  of  the  weight,  Var  (W(n)],  is  needed  to  compute  the  output  power.  The 
weight  variance  is  computed  as  follows.  The  weight  update  equation  is  rewritten  as 
as 


W(n+1)  = W(n)  [\-p|x(n)|2j  + Kd(n)x*(n) 

The  solution  to  this  difference  equation  is  given  by 

n n n 

W(n+1)  = W(o)  n (l -p|x(k)|2)  + p S d(m)x*(m)  11  [l-^x  (k)!2] 

i _ k=m+l  L 
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The  Initial  weight  W(o)  can  be  set  to  zero 
Let 


u 

- F(m)  = d(m)x*(m)  J~|  [\-n)x(k)|2J 


k=m+l 


Averaging  to  obtain  the  mean  weight  yields 


u II  ll 

E[W(n+l)]  = h ^ E[F(m)]  = E[d(m)x*(m)]  ["]  E J^l-^jx  (k>|  2J 
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where  it  is  assumed  that 


E[d(m)x  (m)]  = * * 

E jjx  (m)|  2J  = ,22 

For  d(m)  = ol(m)  + n^m) 

and  x(m)  = I(m)  + n2(m) 

where  and  n2  are  independent, 
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2 2 2 
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The  summation  can  be  expressed  in  closed  form  as 
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In  steady  state, 


E[W(n+l)]  - -2^-2 


<r.”  + < r 
I n 


To  obtain  the  variance  the  mean  squared  weight  is  needed. 
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The  single  sum  is  treated  first. 
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p=m+l  q=m+l 
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p=m+l 

The  expectation  separates  because  the  terms  in  the  product  involve  inputs  at  later  times, 
which  are  uncorrelated.  Using  the  gaussian  assumption, 
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The  double  sum  term  in  E [w^  (n+1)]  involves  E[F(m)F(q)]  which  is  treated  as  follows 


( * * 11  n ) 
E[F(m)F*(q)]  = E jd(m)x*(m)d(q)x*(q)  j”|  \\  (p) t 2 J [l-|Jx  (v)|2J* 


p=m+l  v=q+l 
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In  the  double  sum,  q<m,  so  that 


n n n n m 

n n < > - n n n < »• 

p=m+l  v=q+I  p=m+l  v=m~l  v=q+l 
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Thus,  since  both  p and  v range  over  indices  that  are  greater  than  m and  q, 
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d(m)x*(m)d(q)x*(q)  ["[  l-H-|x(v)|2j 


v=q+l 


The  first  expectation  is  identical  to  e[f2  (m)],  evaluated  previously.  The  second 
expectation,  since  q is  always  less  than  m,  is 
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The  expectation  in  the  center  is 

1 * 3 ) 2 2/  2 2 \ 
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Combining  the  above, 
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The  double  sum  term  is  then 
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Thus  in  steady-state. 
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To  obtain  the  variance,  subtract  the  mean  steady-state  weight  scared,  with  the  result 
that 
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For  a special  case,  if  o-j  »<rQ 
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To  determine  the  cancellation  level,  the  above  result  is  used  in  the  evaluation  of  the 
mean  squared  error,  i.e. , the  power  in  the  canceller  output. 
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If,  In  addition,  \x<r*  «.  1,  then 


E f( « (®)|2  * (1  +|*|2>  <T 

* to  m 


and  the  common  waveform  has  been  cancelled  down  to  the  uncorrelated  noise  flow. 


For  the  tracking  problem,  the  low  signal-to- noise  ratio  case  is  more  critical. 
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The  time  response  la  given  by  the  time  constant,  t , for  the  LMS  algorithm  (1) 

c 


which  can  be  approximated  by 


T * -A_ 

0 2pPin 


where 


A * time  between  iterations 


Pin  * total  input  power. 

These  results  can  now  be  used  to  evaluate  tracker  performance  in  Appendix  IV. 
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APPENDIX  IV 


VARIANCE  OF  THE  SPLIT  BEAM  ADAPTIVE  BEARING  ESTIMATOR 
AND  COMPARISON  TOCRLB 


With  reference  to  Figure  6,  for  the  continuous  adaptive  filter,  It  Is  assumed 

that 

XL  (t)  ■ a(t  - yilne) +n  L (t) 

XR  (t)  =-  sft)  + nR(t) 

where  s(t),  nL(t),  and  nR  (t)  are  zero  mean,  white  random  prooesses  on  [u>L,  with 
power  Pg,  PN,  and  PN,  respectively. 

Also 

d * distance  between  phase  centers  of  half  arrays 
c * speed  of  sound 

6 ■ angle  between  plane  wavefront  and  axis  of  phase  centers 

The  length  of  the  adaptive  filter  is  assumed  to  be  long  In  comparison  to  the  correlation 
times  of  the  random  processes  involved. 

Let  the  steady-state  Impulse  response  of  the  filter  be  h(t)  and  the  transfer  function 
Htco).  The  mean  value  of  the  Impulse  response  Is  the  Wiener  filter,  so,  from 
Appendix  HI,  the  corresponding  mean  transfer  function  Is 

P -Ju^sine 

E(H(«)1  « p-  f -p-  e 0 
*S  *N 

and  the  correlation  of  the  transfer  function  Is 
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The  eagle,  0 , Is  to  be  estimated  by  determining  the  value  of  r for  which  the  Impulse 
response,  h(r),  is  a maximum.  This  is  an  estimate  of  the  delay  between  arrivals  of 
the  signal  wavefront  at  the  half  array  phase  o enters,  and  can  be  converted  to  a bear- 
ing estlmage,  6 , using 


r ■ £ Sin  e . 
c 

The  value  of  § for  which  h(r)  is  peaked  corresponds  to  the  value  for  which  the  derivative 
of  h(r)  with  respect  to  0 Is  sero.  In  the  neighborhood  of  the  tero  crossing  the  derivative 
of  the  Impulse  response  is  approximately  linear.  The  fluctuations  in  the  derivative  of 
the  transfer  function  then  map  through  the  linear  frnotion  to  provide  an  estimate  of  the 
mean  square  error  in  the  angle  estimate.  Thus  the  mean  value  of  the  derivative  of 
the  transfer  function,  the  standard  deviation  of  the  derivative  of  the  transfer  frnctlon, 
and  the  slope  of  the  mean  of  the  derivative  of  the  transfer  function  are  recpired,  at 
the  point  where  6 - 0 . Then  the  errors  in  the  transfer  function  can  be  mapped  Into 
the  errors  in  one' s ability  to  extract  the  peak  of  the  transfer  fraction  as  follows: 


Varl/2  [§1 


L J 


It  --  sin 0 
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This  approach  is  motivated  by  the  treatment  of  BDI  spilt  beam  trackers  by  MacDonald 
and  Schulthelss  In  Reference  [1].  The  derivative  of  the  Impulse  response  and  its 
statistics  are 
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where  the  bend  Is  assumed  to  be  limited  to  (wl.wu>*  The  mean  value  of  the  derivative 
passes  through  sero  when  r»  d/c  sin  6 , so  that  determination  of  the  zero  crossing  pro- 
vides an  unbiased  estimate  of  6 . The  variance  of  the  derivative  Is  needed  as  well. 


Note  that  the  mean  value  of  the  derivative  is  zero' when  e * 0 so  that  the  variance  is 
the  mean  square  value  at  this  point. 


where  v » r - d/c  sin  0 . Using  the  correlation  of  the  transfer  function,  the  variance  is 


where 


The  second  integral  is  zero,  when  v ■ 0,  and  the  delta  function  reduces  the  first 
double  integral  to  a single  Integral,  with  the  result  that 


where  again  the  processes  are  band  limited  to  (oj^,  u^).  Since  T*  d/c  sin  e , 
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The  denominator  of  the  expression  for  Var  [0]  is 
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Thus, 
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The  behavior  with  extreme  cases  of  signal-to-noise  ratio  is  of  Interest.  For  Pg  » PN, 
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and  for  Pg  « P^ 
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The  critical  concern  ia  the  behavior  of  the  variance  at  low  aignal-to- noise  ratios. 

The  objective  therefore  is  to  compare  the  above  expression  with  a performance 
measure  which  is  invariant  with  processor  structure.  The  Cramer-Rao  Lower 
Bound  (CRLB)  is  such  a measure.  The  last  equation  in  this  section  is  therefore  related 
to  .the  CRLB  in  the  next  section. 


COMPARISON  WITH  THE  CRLB 


The  Cramer-Rao  Lower  Bound  (CRLB)  provides  a lower  bound  on  the  variance 
attainable  by  any  unbiased  estimate  of  a parameter.  It  is  computed  in  Reference  (l] 
and  repeated  here  for  comparison  with  the  variance  obtained  by  the  adaptive  filter 
bearing  estimator.  From  equation  (19)  of  [l]. 


Var  1/2  [ 6 ] 2 
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For  small  signal-to-noise  ratios,  i.e.,  Pg  « PN, 
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where  a total  observation  time  of  T seconds  is  assumed,  Pg/P^  is  the  signal-to-noise 
ratio  at  the  split  array  output,  and  d is  the  distance  between  half  array  phase  centers. 

Since  uP^  applies  to  the  discrete  version  of  the  adaptive  filter  and  the  CRLB  uses 
continuous  data  we  must  relate  the  discrete  adaptive  filter  weight  updates  to  a continu- 
ous process.  Let  and  correspond  to  the  feedback  coefficients  in  the  discrete 
and  continuous  versions  respectively.  In  discrete  form 

W((n  + 1)A)  * W(nA)  + < (nA)  x (nA) 
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Taking  appropriate  limits  provides  the  differential  equation  form,  i.e. , 


Thus  for  the  continuous  case  3 Md/A.  This  can  also  be  interpreted  as  though  the 
power  spectrum  in  the  continuous  case  is  the  discrete  power  spectrum  times  A,  with 
H fixed  in  both. 

One  other  term  must  be  reconciled  for  the  two  expressions  to  be  compared,  the 
relation  between  and  the  value  of  E[II(w)]  at  t = t . Previously  it  was  assumed 
that  the  mean  transfer  function  E[H(u)]  had  reached  a steady  state  value  when  deter- 

l/o  ^ 

mining  Var  [0]  above.  Actually  ECH(w)]  is  a function  of  t and  the  expression  for 

l/O  ^ 

E[H(w)]  used  in  computing  Var  [0]  above  is  E[H  (ju,  ® )].  We  desire  ECH(w,  t)] 
where  t « 0 is  when  the  filter  is  turned  on.  It  is  easy  to  show  that  the  continuous  time 
approximation  to  E[H(u>,  t)]  is 


where  T is  the  time  constant  associated  with  the  mean  weights, 
c 

Furthermore  we  assume  that  the  variance  of  the  weights  has  already  reached 
steady-state  so  that  no  correction  for  a transient  variance  term  is  required.  From 
Appendix  m,  it  can  be  seen  that  the  time  constant  of  the  variance  is  approximately 
half  as  large  as  the  time  constant  of  the  mean  for  n ^2  + o-^2)  <<  1. 


Thus,  taking  into  account  the  transient  mean  weight,  for  the  continuous  comparison, 
the  variance  of  the  adaptive  filter  bearing  estimate  becomes 


The  convergence  time  for  the  adaptive  process,  i. e. , the  time  for  the  root-mean- 
square— error  learning  curve  to  decay  by  l/e  of  its  initial  value.  Is 

T - A 

T0  -3^ 

so  that 

-1/2 

v«1/2[«i -Is 

2 d cos  e *S 

c 


r 

CO  »-5 
3 

1 

3 

r -t/Tc  i 

1 . A L 

3 J 

1 • c 

a - 

Comparison  with  the  CRLB  at  low  signal-to- noise  ratios  results  in  the  following  ratio 


The  optimum  selection  of  Tc  involves  the  trade-off  between  weight  vector  noise  and 
mean  value.  Table  10  shows  the  ratio  for  various  T/Tc.  The  minimum  is  at  approxi- 
mately T/Tc  = 1 and  is  1.13.  Hence,  the  split-beam  adaptive  bearing  tracker  is  only 
about  0. 5 dB  degraded  from  the  CRLB. 

It  can  thus  be  concluded  that  the  adaptive  filter  can  be  used  to  estimate  bearing  in 
close  proximity  to  the  theoretical  limits  in  performance.  Although  its  performance 
degrades  at  low  signal- to-noise  ratios,  so  does  that  of  all  processors.  Indeed,  the 
split  beam  adaptive  filter  tracker  behaves  just  as  the  CRLB,  throughout  the  entire 
range  of  signal- to-noise  ratios  and  bandwldths. 
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APPENDIX  V 

DELAY  ESTIMATION  WITH  A DISCRETE  ADAPTIVE  TRACKER 


I 


I 


In  the  preceding  section,  an  expression  for  the  variance  of  the  bearing  estimate  of  a 
continuous  adaptive  filter  configured  as  an  adaotive  tracker  was  developed.  The  esti- 
mate is  based  upon  determination  of  the  peak  of  the  continuous  impulse  response  of  the 
adaptive  filter.  In  practice,  however,  the  adaptive  filter  is  discrete  in  time  and  of 
finite  length,  and  the  peak  of  the  impulse  response  must  be  determined  by  interpolating 
between  the  discrete  sample  points.  The  interpolation  process  increases  the  variance 
of  the  bearing  estimate  with  respect  to  that  of  the  continuous  case.  This  section  de- 
scribes a numerical  method  of  determining  the  peak  of  an  interpolated  Impulse  re- 
sponse, and  develops  an  expression  for  the  variance  of  the  estimate  using  this  method. 

Let  the  discrete  adaptive  filter  weights  be  denoted  h(mT  ) for  m - 0,  1,  ...  M-l,  with 

s 

T the  sample  interval  and  M the  number  of  taps  in  the  time  domain  adaptive  filter. 

The  value  of  the  weight  vector  at  some  time,  t,  between  filter  taps  be  interpolated 
as 

sin  2 B(t-mT  ) 

hj(t)  = 2 h<mT8)  ' 2 B(t-mT  ) (1) 

m=0  8 

This  is  a truncated  expansion  of  the  discrete  samples  in  series  of  orthonormal  func- 
tions. If  the  h(mT  ) are  samples  of  a perfectly  bandlimited  impulse  response,  h(t), 
s 

then  (1)  can  be  considered  as  a finite  approximation  to  the  sampling  theorem. 

In  the  adaptive  tracking  system,  we  are  only  Interested  in  the  location  of  the  peak  of 
hj(t),  not  in  the  entire  function,  so  a numerical  method  is  used  to  locate  the  peak. 

The  peak  of  hj(t)  will  be  in  the  vicinity  of  the  h(mTg)  with  the  largest  value,  so  if  the 
largest  weight  occurs  for  m = q,  then  the  initial  guess  for  the  location  of  the  peak  is 

l0  = «Ts 

The  slope  of  hj(t)  will  be  zero  at  the  maximum,  so  a bisection  technique  is  used  to 
find  the  zero  crossing  of 

V-l 


M-l 

h'j(t)  - £ h(mTs) 
m*0 


cos  2 B(t-mT  ) sin  2 B(t-mT  ) 

B 


t - mT. 


2 B(t-mT  ) 

B 


(2) 


By  observing  the  sign  of  h'j(t),  at  the  current  estimate,  tj,  a new  estimate  is 
computed  as 


■ 'i*,[hw]*t 


i-t-i 
where 

5^  = current  step  size 


(3) 


and 


[hw] 


r+l  if  h'j^)  > 0 


■1  if  h'j^)  < 0 


Each  time  hyt^)  and  b’^tp  have  opposite  sign,  the  zero  crossing  has  been  passed, 
and  the  step  size  is  divided  by  some  F > 1.  The  procedure  stops  when 


-«<  h'jtt^s* 


(4) 


in  which  case  t^  is  the  estimated  location  of  the  zero  crossing,  and  hence,  the  esti- 
mated delay  between  the  tracker  inputs. 


The  variance  of  this  estimate  can  be  determined  in  the  same  way  as  that  using  the  con- 
tinuous adaptive  tracker  in  the  preceding  section.  When  the  tracker  is  implemented  in 
the  frequency  domain,  as  shown  In  Figure  5,  the  mean  and  variance  of  the  complex 
weight  In  the  kth  frequency  bin  has  been  shown,  neglecting  windowing  effects, 
which  are  unique  to  the  frequency  domain  Implementation,  to  be 


E[*k]  - ps  *p„ 


-iff* 


e 


(5) 


and 


Var 


1/2 


M 


MMPn(2Ps+Pn) 

(Ps+Pn)(2-2MM(Ps+Pn)) 


1/2 


(6) 
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■MMMNMI 


■Mi 


We  let  H(k)  = Wfc,  so  that  the  time  domain  weights  can  be  written  as 
M-l 

h<mTs)  = H Z H(k)  e 

k=0 


(7) 


Now,  consider  an  Interpolated  impulse  response,  hj(t),  using  a general  interpolation 
function,  f(t),  that  is 

M-l 

hr(t)  = h(mTs)  f(t-mTs) 
m=0 


M-l 

il 

m=0 


M-l 

I 

k=0 


H(k)  f(t-mT  ) e 
s 


M 


(8) 


Following  the  method  described  for  the  variance  of  the  estimate  with  the  continuous 
adaptive  filter,  the  variance  in  the  determination  of  the  peak  of  hj(t)  can  be  computed 
as 


Var 


1/2 


Var1/2 

' dhj(t)  1 

dt 

dE 

dhj(t)  ' 

dt 

’ J 

dt 

(9) 


Now,  assume  that  the  inputs  to  the  adaptive  filter  are  bandlimited  to  some  bandwidth, 
Bg,  so  that 

H(k)  = 0 for  k > (J-l)  = jBgTgj  ^ (10) 

where  jxj  is  the  largest  integer  less  than  or  equal  to  X.  Further,  in  order  for  h(mTg) 
to  be  real,  we  require 

H(k)  = H*(M-k)  (11) 

Using  (8),  it  can  readily  be  verified  that 
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M-l 


M-l  M-l 


Var 


,1/2 


[ h,i(t)  ] - 3 2 v,r  CH<k,J  2 2 

M k=0  nyO  m2=0 


f’(t-n^Ts)  e 


iT<mrm2>k 


(12) 


and 


dE 


bt 


M-l  M-l 

= ^ S E CH(k)]  £ e 

k=0  m=0 


jfmk 


(13) 


We  want  the  mean  and  variance  of  the  weights  in  (5)  and  (6)  to  meet  the  constraints  of 
(10)  and  (11).  This  leads  to 


P -1— Ak 

*S  JMaK 
e 


PS+PN 


E[H(k)]  = < 


, 0 < k < J-l 
f J— 1 < k ^ M-J 


(14) 


PS  -j^(M-k) 

e M , M-J  < k < M-l 


V PS+PN 


and 


Var  [H(k)]  = 


MMPN(2PS+PN) 


(Ps+Pn)(2-2MM(Ps+Pn)) 


, 0 < k < J-l  or  M-J  s k £ M-l 


(15) 


V 0 

Using  this  in  (12)  and  (13),  we  obtain 


, J-l  < k < M-J 
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'■  1 ■■  ■■'"mi  m*m 


>*PN  (2PS'fPN)(PS't'PN) 
p 2 2-2MM(Ps+Pn) 

S 


1/2 


M-l  M-l 


1/2 


M 


I 

m^O 


I r*r"iV  ‘W.’ 


m2=° 


Jf 

\K=0 


M-l 


m=0 


4f 

\K=0 


j|f(m-A)k 


1/2 


1 

vs 


MPN  (2ps>pN)(ps-pN) 
p 2 2-2/iM(Pg+PN) 

s 


*1 


(16) 


Note  that  depends  upon  the  input  bandwidth,  and  the  choice  of  interpolator  function, 
f(t),  but  not  on  the  signal  to  noise  ratio.  This  is  analogous  to  the  term 


-1/2 


3 „3 

(!)  -4). 

u L 


in  the  expression  for  the  variance  using  a continuous  adaptive  filter.  Equation  (16) 
can  be  further  simplified  by  writing  the  summations  on  k in  closed  form  as 


Re 


J—  1 m ki. 

ejM(mrm2)k 


K=0 


sin 


[: mF  (mrm2>  ] 

[i,(mrm2)] 


cos 


[(iff)  * <mrm2»  ] °7) 


and 


Re 


j|f<A-m)k 


sin 


(18) 
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M_l  (sln 

£ f,,(tz-“T,)< cos  \{^u)  ” (A-m)l 

m=0  I sin  -j^(A-m) 


For  the  interpolation  function  used  here,  <1),  we  have 
cos  2irB(t  -mT  ) sin  2»rB(t  -mT  ) 


f’(t  -mT  ) * — 

1 * ‘s'  (tfmTs) 


z s' 
\ 2 


2jrB(tz-mTs) 


f'(0)  = 0 


f"(tz-mT8) 


2 sin  2irB(t  -mT  ) 2 cos  27rB(t  -mT  ) 


<V"V 


2»B(t2-mT8f 
2rB  sin  2nB(t  -mT  ) 


f”<o>  * - 

Using  (19),  (20),  and  (21),  Kj  can  be  evaluated  numerically  on  the  computer  for  any 
choice  of  signal  bandwidth,  Bg,  interpolator  bandwidth,  B,  filter  size,  M,  and  sample 
Interval,  Tg.  The  computed  value  of  Kj  can  then  be  used  to  compute  the  variance  of 
the  delay  estimate,  tz,  from  (16),  for  selected  slgnal-to-nolse  ratio  and  /*.  Note  that 
for  typical  values  of  /i,  M,  and  (Pg+Pj^).  the  factor 


1 mPN  <2Vpn><Vpn> 

VI  P 2 2-2MM(Ps+Pn) 
s 


in  (16)  is  essentially  independent  of  M,  the  filter  length. 
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It  must  be  noted  here  that  this  result  applies  only  when  the  variance  of  hT(t)  falls  in  the 
linear  portion  of  h'j(t)  near  the  peak  of  E[hj(t)].  This  is  equivalent  to  saying  that  the 
largest  of  the  discrete  weights  must  be  In  the  vicinity'  of  the  true  peak.  This  will  be 
true  when  the  mean  of  the  largest  weight  is  large  in  comparison  to  the  variance  of  the 
weights.  For  a particular  input  signal  to  noise  ratio,  this  can  be  assured  by  selecting 
the  feedback  coefficient,  ft,  sufficiently  small. 

A number  of  simulations  of  an  adaptive  tracker  using  this  interpolation  scheme  have 
been  performed  using  the  configuration  shown  in  Figure  13.  Using  the  results  In  this 
section,  the  estimation  performance  has  been  calculated  on  the  computer,  and  is 
shown  in  the  graph  of  Figure  11  as  a function  of  SNR  and  for  several  values  of  fx.  Also 
shown  are  sample  variances  of  the  estimates  derived  from  the  computer  simulations 
described  in  section  6.  It  can  be  seen  that  when  \i  is  sufficiently  small,  the  theo- 
retical predictions  and  simulation  results  agree  quite  well.  However,  when  fi  is  not 
small  enough  to  be  consistent  with  the  assumptions  of  this  section,  the  variance  of  the 
estimate  is  much  larger.  This  is  because  fluctuations  in  the  weight  vector  cause  the 
largest  weight  to  appear  in  a bin  other  than  the  peak  of  the  mean  weight  vector. 


In  order  to  select  fi  sufficiently  small  to  achieve  the  predicted  results,  we  want  /x  such 
that  the  mean  of  the  largest  weight  is  large  in  comparison  to  the  weight  variance.  Fig- 
ure 12  shows  the  ratio 

E |W 
l_  max_ 

R V2 

Var  W 

as  a function  of  SNR  for  several  values  of  fi,  and  for  the  configuration  simulated  above. 
Comparison  with  results  of  the  simulations  indicates  that  R must  be  nominally  1.5 
to  give  the  predicted  results.  Then,  for  example,  if  performance  were  required  down 
to  -20  dB  for  the  simulation,  fi  < 2~17  would  be  required. 


In  order  to  give  some  feel  for  the  meaning  of  these  results  in  terms  of  bearing  estima- 
tion, Figure  11  also  shows  predicted  variance  of  a bearing  estimate  for  a broadside 
target  for  two  arrays,  a small  tactical  array  with  11.46  feet  between  split  array  phase 
centers,  and  a towed  array  with  114.6  feet  between  phase  centers. 
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APPENDIX  VI 

TRACKING  BEHAVIOR  OF  THE  MEAN  WEIGHTS  FOR  VARIOUS  SIGNAL 
MODELS ^AND  A LINEARLY-TIME-VARYING  DELAY 

1.  Tracking  Behavior  of  an  Adaptive  Filter  TracKer  for  Single  Frequency  Inputs 

with  Linearly  Time-Varying  Bearing  in  Uncorrelated  Noise 
• 

An  LMS  adaptive  filter  is  configured  as  a canceller  between  two  half  array  beamformed 
outputs  to  perform  split  beam  bearing  tracking.  The  tracking  performance  is  analyzed 
for  a narrowband  signal  that  is  moving  such  that  the  delay  between  split  array  phase 
centers  is  linearly  changing  with  time.  The  time  varying  mean  weights  are  derived  in 
both  transient  and  steady  state  conditions  and  compared  with  previous  results. 

The  LMS  adaptive  filter  is  configured  as  a bearing  tracker  in  Figure  2.  The  objective 
is  for  the  adaptive  algorithm  to  properly  estimate  and  track  changes  in  the  propagation 
time  between  split  arrays  and  map  that  time  estimate  into  an  estimate  of  signal  energy 
arrival  angle.  The  narrowband  signal  case  in  uncorrelated  noise  is  treated  herein.  It 
is  assumed  that  0(t)  is  such  that  the  delay,  r(t),  between  array  phase  centers  is  linear 
with  time.  Since 

r(t)  = sin  ©(t) 
o 

for  a split  array  aperture  of  d meters  in  a sound  speed  of  C©  meters/sec,  the  assumption 
of  linear  r(t)  corresponds  to  assuming  small  linear  ©(t),  such  that  sin  ©(t)  » ©(t)  and 

-(t)  a 7T-  ©(t)  - ct. 

°o 

1 
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The  model  for  the  moving  case  is  as  follows.  The  output  of  the  upper  beamformer  in 
Figure  2 produces  the  reference  waveform,  d(t),  where 


d(t)  = A cos  (w0tf  <p)  + ^(t) 

The  filter  input  is  the  lower  beamformer  output 
x(t)  = A cos  (wQt  - r(t)  + <P)  + n2(t) 


* A cos|u)0t  - -g—  sin  6(t)  + <f>j  + n2(t) 


= A cos  ; t/(t)  + n2(t) 

where  if(t)  is  the  total  signal  phase.  The  frequency  of  the  signal  portion  of  x(t)  is  the 
time  derivative  of  so  that 


2*f(t)  = = w0  fl  - ^ cos  ©(t) 


Thus  the  effect  of  the  moving  plane  wave  signal  is  to  introduce  a frequency  shift  on  the 
sine  wave.  The  ability  of  the  adaptive  filter  to  track  this  moving  frequency  shift  will 
now  be  analyzed. 


The  algorithm  for  changing  the  complex  weights  in  the  adaptive  filter  is  given  by 
W(n+1)  = W(n)  * (d(n)  - XT(n)W(n)]  X*  (n) 


where 


= W(n)  + m (d(n)  X*(n)  - X*  (n)XT(n)W(n)] 


= filter  weight  vector  a time  sample  n 
= desired  signal  at  time  sample  n 


observed  data  vector  of  samples  within  the  tapped  delay  line 


It 


and  where  the  symbols  * and  T denote  complex  conjugate  and  vector  transpose  respec- 
tively. Figure  8 shows  the  data  vector  X(n)  in  the  tapped  delay  line.  The  scalar  d(n) 


ju  nAt 

d(n)  = tfg  e + n (nAt) 


and  the  vector  X(n)  is 


jw  (nAt  - r(nAt)) 


j uQ  (nAt -6  - r(nAt-6)) 


(nAt-26  - r(nAt-2<5) 


X(n)  = cr 


ju.Q  (nAt-M<5  - r(nAt-M<5)) 


where 


6 = time  delay  between  taps 
At  = algorithm  sampling  time  (usually  At  = 6) 
M = number  of  taps. 


The  equation  of  the  mean  weight  vector  is  obtained  by  averaging  (1)  and  assuming  that 
the  data  sequence  X(n)  is  independent  over  time. 

E[W(n+l))  = E(W(n)]  + n I R^n)  - R^n)  E[W(n))  | (4 
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Using  the  inputs  in  (2)  and  (3),  R^  and  R^  are  calculated  as  follows: 


Rdx(n)  * E[d(n)  X (n)]  = a S 


J«0  r(nAt) 


[r(nAt-6>+(5] 


jc^  [-r(nAt-M6)  + M<5] 


(5) 


R (n) 
xx  ’ 


, S - r(nAt-16)  + (k-l)6] ) 

« E[X(n)X  T(n)J  - cr*  e 0 J 


+ a 2 I 
n 


(6) 


where  k and  1 range  from  O to  M and  I is  the  MxM  identity  matrix  for  the  covariance  of 
the  uncorrelated  noise. 

Note  that  both  second  order  statistics  are  dependent  on  nAt,  reflecting  the  non-station- 
arity  of  the  input  processes.  Define  the  matrix  D (n)  as  follows: 

I Ju  r(nAt) 
e 


D(n) 


ju;o  [r(nAt-<5)  + 6] 


JuiQ  [T(nAt-M<5)  + M6] 


(7) 
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Then 


Rdx(n)  “ & SD(n) 


R^n)  »<r2s  D(n)  D^n)*  I 


(8) 

0) 


If  -(t)  * ct,  as  assumed  at  the  outset,  then 


\u  C oAt 
/ o 


D(n)  = 


J*0  (l-c)<5 


jo;  (1-c)  M6 


and 


JUo  CnAt  D(0) 


(10) 


Rdx(n)  sffSe 


2 J^0  c nAt 


D(0) 


Rxx(n)  = (Tg  D(0)  D(0)  *T  + ct2  I 


(U) 

(12) 


Using  (11)  and  (12)  in  (4),  the  mean  weights  are 


E(W(n+l)]  = J I - M [an2  I +a2g  D(0)D  T(ojJ  J E(W(n)J  «■ 


2 iw0  c nAt 


D(0) 


(13) 


The  difference  equation  (13)  can  be  solved  explicity  for  E[W(n)].  Assuming  that 
E(W(0)]  * 0,  i.  e.  zero  initial  conditions, 


9 jw  c (nf  l)At  . 

E(W(n)l  = u a g e ° [I-A]"a  (I-A11]  D(0) 


(14) 


where 


A • e 


-lu.  cAt 
J o 


[i-  u(cn2  I+<r2gD(0)D(0)*T)] 


(15) 
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then 


E[W (n)]  = a_(n+1)  ^ 


jwQ  (l-c)«5 


jwo  (l-c)M<5 


(22) 


The  weights  are  all  of  magnitude 


3 

i>MT 


2 

V a sa 


1 - a+an  ^aQ^+Mffg  ^ 


(23) 


but  have  a phase  that  is  precessing  across  the  tapped  delay  line  as  time  progresses.  The 
weights  may  be  viewed  as  a sinusoid  at  the  shifted  frequency  of  the  input,  uq  (1-c),  with 
a phase  that  is  "barber  polling"  along  the  delay  line.  The  weights  are  shown  in  Figure  9 
as  time  changes. 

2 2 

In  the  static  case,  c = 0,  or  = 1,  3 = a e/cr  , and 

O U 

/ ' 


E[W(~)J 


2 , 2 

1+M  4 \2 


jo;  M6 


(24) 


which  agrees  precisely  with  previously  derived  results  for  the  LMS  algorithm  in  steady 
state  with  a staionary  sine  wave  input  in  uncorrelated  noise. 
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Conclusions 

The  transient  and  steady  state  tracking  behavior  of  the  mean  value  of  the  LMS  adaptive 
filter  tracker  weights  has  been  derived.  The  filter  corrects  a frequency  shift  and 
doppler  precession  due  to  the  angular  movement  of  the  target  in  bearing. 

2.  Tracking  Behavior  of  an  Adaptive  Filter  Tracker  for  a Broadband  White  Signal 
with  Linearly  Time-Varying  Bearing  in  Uncorrelated  Noise 

The  adaptive  filter  split-beam  bearing  tracker  analysis  is  now  extended  to  include 
a broadband  spectrally  white  input  signal  that  is  moving  in  bearing.  The  time  varying 
mean  weights  are  derived  in  both  transient  and  steady  state  conditions. 

The  LMS  adaptive  filter  is  configured  as  a bearing  tracker  in  Figure  8.  The  objective 
is  for  the  adaptive  algorithm  to  properly  estimate  and  track  changes  in  the  propagation  time 
between  split  arrays  and  map  that  time  estimate  into  an  estimate  of  signal  energy  arrival 
angle,  0(t).  It  is  assumed  that  0(t)  is  such  that  the  delay,  r(t),  between  array  phase 
centers  is  linear  with  time.  Since 
d 

r(t)  = — sin  0(t) 
uo 

for  a split  array  aperture  of  d meters  in  a sound  speed  of  C0  meters/sec,  the  assumption 
of  linear  r(t) corresponds  to  assuming  small  linear  0(t),  such  that  sin  0(t)  » 9(t)  and 

r(t)  a 9(t)  = ct. 
o 

The  input  signal,  s(t),  in  this  case  is  a zero  mean  Gaussian  random  process  that  is 
spectrally  white  over  the  band  corresponding  to  the  sampling  frequency.  The  total  signal 

2 2 

power  is  <r  and  the  total  uncorrelated  noise  power  is  a . As  in  the  prior  analysis,  one 
s n 

split  array  output,  d(n),  is  the  desired  signal  for  the  adaptive  process,  and  the  other  is 
the  input,  x(n),  to  the  adaptive  filter.  The  ability  of  the  adaptive  filter  to  track  the  bear- 
ing changes  in  the  broadband  plane  wave  signal  process  is  analyzed  in  the  next  section. 

Analysis 

The  algorithm  for  changing  the  complex  weights  in  the  adaptive  filter  is  given  by 

W(n+1)  = W(n)  + m [d(n)  - XT  (n)W(n)]  X*  (n) 

= W(n)  + n [d(n)  X * (n)  - X*  (n)XT  (n)W(n)]  (1) 
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where 

W(n)  = filter  weight  vector  at  time  sample  n 
d(n)  = desired  signal  at  time  sample  n 

X(n)  = observed  data  vector  of  samples  within  the  tapped  delay  line 

and  where  the  symbols  * and  T denote  complex  conjugate  and  vector  transpose  respec- 
tively. Figure  11  shows  the  data  vector  X(n)  In  the  tapped  line.  The  scalar  d(n)  is 


d(n)  s(nit  - r(nAt))  + nj  (nAt) 


and  the  vector  X(n)  is 


X(n)  = 


/ s(nAt) 

[ s(nAt  - 6) 
s(nAt  - 26) 


+ N,  (nAt) 


(2) 


where 


6 

At 

M 


^ s(n  At  - M6)  I 

- time  delay  between  taps 
= algorithm  sampling  time  (usually  At  = 6) 
= numbar  of  taps. 


It  is  assumed  that  T(t)  - ct  and  that  the  signal  is  white,  so  that 


E[s(n)s(m)j  = <j“  A (n-m) 

where 

A (n-m)  = I1  lfn=m 

| 0 otherwise 


(3) 


(4) 


1 


i 


The  difference  equation  for  the  mean  of  the  adaptive  filter  weight  is : 

E [W(n+l)J  = [I  - n Retail  E IW(n))  + pR^  (n) 

where 

p • feedback  coefficient 
I » MxM  identity  matrix 
Rdx(n)  = E[d(n)  X*(n)) 

^xx(n)  ~ E[X(n)  X*^(n)J 

If  Rxx<n)  is  independent  of  n,  and  if  the  initial  mean  weight  E[W(o)J  * O,  than 
equation  (5)  can  be  re-written  as 


(5) 


E(W(n)j  = n 


(i  - MR*/-*-1  R^  0» 


<»3) 


The  second  order  statistics  can  be  calculated  from  the  input  waveforms,  as  follows: 


RxxW  = <as  + %2)  1 


I 


Rdx<n)  ' 'i 


A (cnAt) 

A (cnAt  - 6) 
A (cnAt  - 26) 


A (cnAt  - M6) 

Using  (6),  the  mean  weight  vector  at  the  nth  iteration  is  given  by 
o X1  „ 2 n-k'1  / A(ckAt) 

L,  [i-^;+oi  / 

k=o 


E[W(n)]  = na‘ 


A(ckA-6) 


(7) 


(8) 


0) 


A(ckAt  - M6) 


In  the  static  case,  c ■ o and  only  the  first  weight  is  non-zero,  with  mean  value 

E (WQ(n)J  = - 3 [ 1 - (1  - #*(cr 2 + o “ ))  J (10) 

cr  + a ° u 

s n 

which  converges  to  the  Wiener  Filter  for  the  broadband  stationary  case  as  n - 

In  the  dynamic  case,  the  weights  are  a moving  set  of  spikes  that  are  changing  with  amplitude 
as  the  signal  moves  and  as  the  weights  converge.  The  total  weight  vector  is  the  sum  of  the  vectors 
in  (9).  The  weights  can  be  viewed  as  a sliding  window  of  exponentially  growing  responses,  or  as 
a moving  weight  at  the  leading  edge  that  leaves  behind  it  an  exponentially  decaying  wake. 

This  can  be  seen  by  examining  the  weight  at  die  leading  edge  of  the  response.  In  (9), 
the  leading  edge  will  occur  at  the  latest  time.  If  the  filter  is  sufficiently  long  so  that  the 
response  still  falls  within  the  tapped  delay  line,  i.  e. , M > n,  then  the  amplitude  and 
location  of  the  leading  weight  are  found  by  examining  the  term  in  the  summation  for  which 
^ = n !•  The  amplitude  of  the  leading  edge  is  per,?  and  its  location  is  at  tap  number 
c(n-l)  (assuming  that  A t = 6).  If  c = 1 then  the  signal  moves  one  tap  per  iteration  and 
the  adaptive  filter  tracks  the  movement.  If  the  signal  changes  more  slowly,  then  c «.  1 
and  the  leading  edge  moves  more  slowly  then  the  iteration  rate. 


For  the  special  case  where  c = 1 the  weight  vector  in  (9)  can  be  readily  expanded, 
as  follows,  letting  r = 1 - p (cr2  + a 2 ^ 


Equation  (12)  shows  the  decaying  wake  behind  the  leading  weight  which  shifts  along 
the  delay  line  as  n increases.  Figure  10  shows  this  effect. 

For  c < 1 the  weights  move  more  slowly  and  the  basic  model  herein  tends  to  become 
less  realistic.  The  signal  sequence  used  is  totally  uncorrelated  in  time.  In  general  this 
is  not  the  case.  The  impact  of  this  assumption  is  to  have  weights  respond  only  at  the 
exact  correct  alignment  of  input  delays  and  tap  delay  values.  In  the  band  limited,  non- 
white signal  case,  correlation  will  exist  even  at  non-integer  delay  shifts,  and  larger 
weight  responses  should  be  expected  at  the  leading  edge  for  slower  moving  signals.  It 

is  shown  in  the  next  section  that  the  amplitude  of  the  leading  edge  decreases  monotonic- 

2 2 2 2 
ally  with  c from  <rg/(<rs+  o^),  (the  value  for  c - o),  to  fiars  , (the  value  for  which 

movement  is  so  fast  that  the  signal  samples  decorrelate  totally  at  each  time  sample). 

The  extent  of  the  wake  and  the  height  of  the  leading  edge  will  depend  on  signal  dynamics. 


Conclusions 

The  transient  and  steady  state  tracking  behavior  of  the  mean  value  of  the  LMS  adap- 
tive filter  tracker  weights  has  been  derived  for  a spectrally  white  process  with  linear 
time  varying  bearing  dynamics.  The  filter  weights  are  a moving  tap  which  leaves  be- 
hind it  an  exponentially  decreasing  wake.  Tracking  the  bearing  involves  estimating  the 
delay  location  of  the  leading  edge  of  the  weights. 

3.  Tracking  Behavior  for  Band-Limited  Broadband  (Correlated)  Signal  with 

Linearly  Time-Varying  Bearing  in  Un-Correlated  Noise  for  Small  Signal-to- 

Noise  Ratios 

The  adaptive  filter  split-beam  bearing  tracker  analysis  is  now  extended  to  Include 
band-limited  broadband  (correlated)  signals  that  are  moving  in  bearing.  The  time- 
varying  mean  weights  are  derived  from  the  transient  and  steady-state  conditions  for 
low  signal-to-noise  ratios. 

The  conditions  of  the  previous  section  hold  with  one  extension.  The  slgn.-vl  spectrum 
is  no  longer  white.  Equations  2 and  3 of  the  previous  section  are  still  valid.  However, 
R^n)  and  R^n)  (Equations  7 and  8 of  the  previous  section)  are  now  given  by 
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r 


Rdx<n>  = 


XX 


where 


\ 


2 2 
% I+(rs  A 


A = 


SYM 


Since  R . (n)  is  independent  of  n,  for[E  W(o)]  = 0 


n-1 


n-k-l 


E W(n) 


= I [i  - Rdx(k) 


k=0 


let 


3 (l  -M^n2)  I- 


A 


and 


D(k)  = 


' p i ck  At  I 
p |ckAt-<5 


pick  At- Mil 
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H 


then 


e[w(b>]  = ti*l  J Rn"k-1  D(k) 

k=0 

2 2 

Assume  small  signal-to-noiae  ratio,  i.e.,  o8  <<on  • 

Then 

Bxx<">  * % 1 

R , (l  - ^ I 

If  we  select  At  = 5,  i.e. , iterate  at  the  sample  rate,  then  the  mean  weight  for  small 
signal-to-noise  ratio  is  approximately 


Assume  that  the  filter  is  long  enough  so  that  M = n-1,  and  let  j Me  so  that  we  are  at 
the  end  of  the  filter.  For  r = 1 - poQ2 

M 

e[wMcI.)]  • rM  | r-*P  c6 

k 0 

M -k 

M<r2rMpM6c  £ (rpc5) 
k=0 
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-M-l 


fKx  2 rM  pMSc  1 - (rP  ) 


1 - <rPc8) 

n*  p*  -M 

= u<r 2 M pMc8  rPCd  - (rPC6) 
s rPc8  - 1 

2 rM  p (M+l)c5  _ x 

rPcA-l 

• 

This  is  the  transient  behavior  of  the  mean  weight  for  small  signal-to-noise  ratios.  In 
steady-state,  i.  e. , M arbitrarily  large,  one  can  readily  see  the  relationships  between 
signal  correlation,  P,  signal  dynamics,  c,  and  algorithm  dynamics,  p. 


lim  E WMc(M) 
M - ® 


2 

1 - rPC 


- [*  - "%] 


There  are  several  cases  of  interest,  Ii  p=  1,  then  the  signal  is  totally  correlated 
from  tap  to  tap,  regardless  of  signal  dynamics.  For  this  case 


E WMc<  > = 7? 

n 


which  is  a small  signal-to-noise  ratio  approximation  to 


2 _ 2 
tr  i-ff 
n s 

which  is  what  one  would  expect  in  the  static  case. 

If  P=  0,  the  signal  is  totally  un-correlated  and  again  signal  dynamics  should  not 
affect  the  result.  For  this  case, 

E(W M„<  )]  = H-O’e 
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which  agrees  with  the  results  of  the  previous  section  (Equation  12)  where  dynamics 
tended  to  de-correlate  the  signal. 

If  there  are  no  signal  dynamics,  i.  e. , c = 0,  then 
ewmo‘">  -^T 

n 

independent  of  P,  which  is  again  the  expected  small  signal-to-noise  ratio  static  result. 

The  tradeoff  between  dynamics  and  signal  correlation  can  be  seen  by  examining 
the  term  Pc°  in  the  expression  for  E[WMc («>)].  As  c decreases  Pc0  looks  more  like 
unity  and  signal  looks  more  correlated.  As  c increases  P (which  is  less  than  one)  is 
raised  to  a higher  power  and  the  signal  has  become  less  correlated. 
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APPENDIX  VII 

Transient  behavior  of  the  LMS  Adaptive  Filter  response  to  variable 

frequency  spectral  lines. 

ABSTRACT 

The  transient  behavior  of  the  LMS  adaptive  filter  is  studied  when  configured  as  a 
canceller  operating  in  the  presence  of  a fixed  or  variable  complex  frequency  sine-wave 
signal  buried  in  white  noise.  For  a fixed  frequency  signal,  the  mean  weights  are 
shown  to  respond  to  signal  more  rapidly  than  to  noise  alone.  For  a chirped  signal,  a 
fixed  parameter  matrix  first-order  difference  equation  is  derived  for  the  mean  weights 
and  a closed-form  steady-state  solution  obtained.  The  transient  response  is  obtained 
as  a function  of  the  eigenvectors  and  eigenvalues  of  the  input  covariance  matrix.  Suf- 
ficient conditions  for  the  stability  of  the  transient  response  are  derived  and  an  upper 
bound  on  the  eigenvalues  obtained.  Finally,  the  mean-square  error  is  evaluated  when 
responding  to  a chirped  signal.  The  gain  coefficient  of  the  LMS  algorithm  is  determined 
that  minimizes  the  mean-square  error  for  chirped  signals  as  a function  of  chirp  rate 
and  signal  and  noise  powers. 


INTRODUCTION 


The  LMS  adaptive  filter  has  been  proposed  and  used  in  situations  where  the  statistics 
of  the  input  processes  are  unknown  or  partially  known  The  structure  of  the  LMS 

algorithm  for  adjusting  the  weights  of  the  adaptive  filter  requires  quadratic  operations 
on  stochastic  input  data  which,  in  general,  are  difficult  to  analyze.  Under  the 
assumption  of  statistically  Independent  data  samples,  the  mean  weight  vector  and  the 
covariance  of  the  weight  fluctuations  have  been  obtained  for  a variety  of  stationary  input 
data  statistics  Special  configurations  of  the  LMS  algorithm,  such  as  noise 

cancelling  ^ line  enhancing^’ 11 1,  spectral  analysis, 12 ^ and  single  fre- 
quency line  detection7-1^,  have  been  studied  in  considerable  detail.  The  special 
characteristics  of  the  LMS  filter  configuration  have  been  used  to  aid  in  the  analysis 
of  the  behavior  of  the  algorithm. 

The  purpose  of  this  paper  is  to  present  some  exact  analytical  results  for  the  LMS 

algorithm  configured  as  an  adaptive  noise  canceller  when  the  input  process  consists  of 

a chirped  sine  wave  in  additive  stationary  <>hite  noise.  .Although  some  previous  work 

f 1 3 • 1 6 ] 

on  LMS  algorithm  behavior  in  a non-stat ionary  environment  has  been  published 
only  one  has  investigated  the  response  of  the  LMS  algorithm  to  chirped  sinusoids 
in  white  noise.  The  analysis  is  performed  by  assuming  the  chirping  is  slow  enough 
so  that  a quasi-stationarv  model  for  the  mean  weights  can  be  used.  In  this  paper, 
exact  analytical  results  are  obtained  for  the  chirped  sinusoidal  signal  with  arbitrary 
chirp  rate.  Since  the  adaptive  cancelling  of  dynamic  signals  is  a key  element  in 
cancelling,  line  enhancing  and  frequency  tracking,  the  analytical  results  for  the  above 
model  have  wide  applicability. 

Two  principal  results  of  this  paper  are 

1.  A closed  forqi  analytical  expression  for  the  LMS  mean  weights  in  a dynamic 
signal  environment. 

2.  Explicit  trade-off  results  between  filter  parameters,  w eight  variances, 
rmean-square-error,  and  input  signal  dynamics. 

*\ 

The  latter  result  is  of  special  interest  since  it  shows  explicitly  the  compromise 
between  fast  adaptation  in  order  to  respond  to  variations  in  the  input  statistics  and 
slow  adaptation  to  reduce  the  fluctuations  in  the  adaptation  process  itself. 
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For  the  narrowband  signal  in  white-noise  case,  the  configuration  shown  in  Figure  40 
can  be  used  to  model  the  above  LMS  algorithm  functions. 

With  reference  to  Figure  40, 

j(w  nAt  + J(nAt)  /2  +0) 
d(n)  = age  0 + n^nAt) 

A is  chosen  so  that  n^nAt)  and  n^nAt  -A)  are  un -correla ted.  On  the  other  hand, 
because  the  desired  signal  is  a chirped  sine-wave,  it  decorrelates  more  slowly  than 
the  noise. 


DYNAMIC  MODEL  FOR  THE  INPUTS 


1 


The  algorithm  for  changing  the  complex  weights  of  the  adaptive  filter  is  given  by  , 

W(n+1)  * W(n)  + h [d(n)  - XT(n)  W(n)]  X*(n) 

* W(n)  + [d(n>  X‘(n>  - X*(n>  XT(n)  W(n)]  (li 

where  W(n)  * filter  weight  vector  at  time  n,  d(n>  * desired  signal,  X(n)  * observed 
data  vector  at  time  n,  and  where  * and  T denote  complex  conjugate  and  vector 
transpose  respectively. 

Averaging  equation  (1)  and  assuming  1)  the  data  sequence  X^nt  is  statistically  inde- 
ll-41 

pendent  over  time  and  2)  the  present  weight  vector  and  the  present  data  vector  are 
statistically  independent^11!  yields 

E(W(n+l)l  = E[W(nil  + »*[Rdx<B>  “ E (2> 

where  R^n)  = E (d(n)X*(n)l,  R^jn)  * E [x*(n>XT  (n)]. 

In  practice  the  algorithm  sampling  interval  (At)  is  usually  chosen  to  correspond  to 
the  delay  *5  between  the  taps  of  the  adaptive  filter.  Furthermore  (At)  is  usually 
chosen  to  correspond  to  independent  samples  of  the  noisy  data.  Hence  the  delay 
is  chosen  to  be  integer  multiples  of  (At)  in  order  for  the  noises  in  the  two  inputs 
to  be  un-correlated.  On  the  other  hand,  the  longer  that  A is  chosen,  the  less 
correlated  is  the  signal  component.  Thus  choice  of  A = 6 is  the  best  that  can  be 
accomplished. 1 

When  the  input  consists  of  a complex  sine-wave  with  linearly-varying  frequency  in 
additive  noise, 


. ? 

t * «Jt  /2  *■  0) 

d(t)  « <7  e ° + n(t)  (3  * 

s 


‘Other  integer  values  of  <5  for  the  bulk  delay  A can  be  studied  using  the  subse- 
quent analysis  and  the  results  show  that  A = a yields  the  best  filter  performance. 
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where  ^s2  ■ signal  power,  uQ  “ signal  frequency,  w * rate  of  change  of  signal 

frequency,  o = random  phase  of  signal  and  the  noise  Is  independent  of  the  signal 
2 

with  noise  power  oQ  and  normalized  covariance  matrix  G,  then 

R (n)  * 2 G + a 2 D(n)  D*(n)T 

xx  n s 


Rdx(n)  = CTs  D(n> 


where 


T - j6  2 2 

D (n>  = e \e  eJ  

JUom6  ejdm<52n  - ejm252I'/2  . . 
e^oM6  eJuM62Q  - ejM2<52w'/2  j 


(4 1 
(5) 


(6) 


with  M = number  of  complex  weights. 


Using  Eqs.  (4)  and  (5)  in  Eq.  (2)  yields 


E(W(n+l)]  = I - p (c  2 G + <7  2 D(n)D*(n)  f)  E(W(n>)  + p cr : 2D(n) 

ns  I S 


For  white  noise,  G = I.  Define 


s T 

M(n)  = I + -=s-  D(n)D*(n) 1 


(T) 


(8' 


For  any  n,  the  eigenvectors  of  M(n)  are  the  vector  D(n)  and  any  set  of'(M-l)  vectors 
orthogonal  to  D(n) . The  associated  eigenvalues  are 
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SOLUTION  OF  EQ  (7)  FOR  THE  MEAN  WEIGHT  BEHAVIOR 

Since  M(n)  is  Hermitian,  there  exists  a unitary  transformation  P(n)  which  diagonalizes 
M(n)  for  each  n, 


P(n)  M(n)  P_1(n)  = X = Diag  (Xj,  X2  • ..Xm', 


The  X.  are  not  functions  of  n.  Due  to  the  special  form  of  D(n), 


T n T 

D1  (n)  = e ° Vn  D (o) 


(10) 


(11) 


where 


' V = Diag  (a 


2 M.  jiJ2 

, a , , a *),  a = eJ 


Also  P(n)  can  be  written  in  terms  of  the  eigenvectors  of  M(n), 
P+(n)  = y^[D(n)'  Rl(n)'  Rm-l(n)] 


(12) 


where+  = conjugate  transpose  and  R^,  Rg  • • • are  ^"1  mutually  ortho-normal 

vectors,  also  orthogonal  to  D(n)  for  each  n.  Using  Eq.  (11), 


P(n)  = P(o)  (V*) 


(13) 


Using  Eqs.  (8),  (11)  and  (13)  and  defining  Z(n)  = P(n)  E[w(n)],  Eq.  (7)  can  be  written 
in  terms  of  Z only  as 


Z(n+1)  = P(0)  V*  P-1(0)  [l-^2x]  Z(n)  + p.cJs2  P(0)  V*  D(0) 


(14) 
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Since  Eq.  (14)  is  a constant  coefficient  linear  difference  equation,  with  P(0)  * P , it 
follows  that 


1 I" 

2(m  - po  V*  P/1  [i  - Z(0) 


E !p.v,v1Ki]j  pc 

m*l  I 1 


V*  D(0) 


Before  investigating  the  general  case  of  Eq.  (15),  consider  the  fixed  frequency 
sinusoid  signal  case  when  V = I and  Eq.  (15)  simplifies  to 


Z(n)  = [i-^o-^]21  Z(0)  +pas2  1 S 


where  ST  = (/mT  0,  0,  ...  0).  Expressing  the  matrix  sum  in  closed  form 


(7  2 -1 

Z(n)  = [l-,on2x]  Z(0)  + -Sg-  X S [i  - (I-Han2  X )n] 

n 


Thus,  using  Eq.  (9),  the  compnents  of  Z(n)  are  given  by 


- 11  - * m‘)]  2i<0)  * Srr  I1  - [ ■ 


2*Mos2j] 


^ l 


J = 2,3 M 
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Hence,  for  ^(on2+  M crg2)  < 1,  the  response  of  the  weights  to  the  signal  frequency  is 
more  rapid  than  to  any  other  frequency.  If  z^(0)  * 0,  j * 1,2, . . . M,  then  z^n)  is  the 
only  response, 


(19) 


Transforming  back  to  the  original  coordinate  system, 


E(W(n)l  = PQ 


Z(n>  = 


(20) 


Hence,' the  mean  weights  are  scaled  versions  of  the  desired  signal  response.  From 

Eq.  (IS),  note  that  the  time  it  takes  the  filter  to  adapt  from  zero  initial  conditions  and 

learn  the  signal  is  less  than  the  time  required  to  forget  the  signal  if  it  disappears. 

That  is,  from  Eq  (19),  if  z.(0)  = 0,  j = 1,2,  ...  M,  signal  response  time  is  propor- 
2 2 * 

tional  to  1 - ^(cr  +Mcr  ).  If  the  signal  suddenly  disappears  so  that  z.(0)  r 0,  then 
n <s 

from  Eq.  (17)  with  cr  4 = 0,  the  decay  time  towards  z.(n)  = 0,  is  proportional  to 

2S 
1 - H0TB*. 
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STEADY-STATE  WEIGHT  BEHAVIOR 


The  explicit  solution  of  Eq.  (14)  requires  evaluation  of  the  eigenvalues  and  eigenvectors 
of  the  matrix  operator  in  brackets  in  Eq.  (15)  (see  Appendix).  However,  the  steady- 
state  solution  to  Eq.  (14)  is  obtainable  without  knowledge  of  the  eigenvalues.  In  Eq.  (15), 
set  Z(0)  = 0 (zero  initial  conditions)  without  loss  of  generality.  Let  Q be  the  matrix  of 


eigenvectors  of  the  matrix  PQV  1 
the  matrLx  of  eigenvalues.  Thus 
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In  Appendix  I,  it  Is  shown  that  | A { |<1,  for  all  i If  0 < n (<rn2+ M crg2)  <2.  Let  ZflS  = 
lim  Z(n).  Then,  using  Eq.  (23),  Eq.  (21)  becomes 

Q ^ oo 

Zoo  = P<r  2 Q(I-.\f 1 Q"1  p0  V*  D(o) 
ss  s y 

^p^v-ii-k^i^dW  (21 

1 + klffg2  Do+  [v  - (W<Tn2)  l]"1  D(0) 


The  steady-state  weights  are  the  quantities  of  interest.  Note  that  they  will  be  time 
varying,  even  though  the  adaptive  filter  Is  in  steacfy-state.  Here,  steady-state  implies 
that  the  adaptive  filter  has  converged,  interpreted  as  the  convergence  of  the  trans- 
formed weight  vector  Z(n).  However,  the  filter  has  converged  to  a time-varying 
solution  to  follow  the  time-varying,  non- stationary  input  signal.  Thus,  when  Zgs  is 
inverse  transformed  back  to  the  mean  value  of  the  weights,  the  transform  is  via  the 
eigenvectors  of  the  input  covariance  matrix,  which  are  time-varying.  Let  E[\Vggtnt] 
denote  the  mean  value  of  the  steady-state  weights  at  time  n. 


e[wm<»>]  - P'1!"'  Z 
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As  a check,  Eq.  (25)  can  be  compared  with  the  steady-state  value  of  the  weights  in 
the  stationary  case,  i.  e. , with  u>  = 0.  For  that  case  V = I,  D(nl  = D(0)  and 


n 


which  agrees  with  Eq.  (20)  when  n — » . 


(26  > 


VII- 11 


Ji 


r inmum 


Computer  evaluation  of  the  steady- state  mean  weights  In  Eq.  (25)  Is  presented  In 
Figures  41-44  for  f = 5 Hz/sec2  and  in  Figures  45-48  for  f * 1.25  Hz/sec2.  In  all  cases, 
the  filter  has  128  taps  with  m **  . 1.  The  signal-to-noise  ratios  are  varied  from 
unity  to  10~2.  The  figures  display  the  magnitude  and  phase  of  the  weights  across  the 
filter.  Three  interesting  phenomena  are  displayed  in  these  figures: 

1.  As  the  signal-to-noise  ratio  decreases,  the  adaptive  filter  uses  more 
of  the  taps  but  at  lower  amplitudes, 

2.  The  tap  phases  follow  the  movement  of  the  linearly  varying  frequency 
input, 

3.  As  f Increases,  the  taps  at  the  far  end  of  the  line  contribute  relatively 
less  to  the  filter  output  than  those  taps  at  the  beginning  of  the  line. 

These  phenomena  can  be  explained  as  follows: 

1.  The  two  sources  of  randomness  that  contribute  to  the  filter  output  mean- 
square-error,  are  input  noise  and  algorithm  noise  (weight  misadjustment).  The 
contribution  of  the  input  noise  to  the  mean  square  error  decreases  linearly 
with  the  number  of  taps  whereas  the  algorithm  noise  increases  linearly 
with  the  number  of  taps.  Thus,  at  high  input  signal-to-noise  ratios,  the 
algorithm  noise  is  the  limiting  factor  and  few  taps  are  needed.  At  low 
input  signal-to-noise  ratios,  input  noise  is  the  limiting  factor  and  a 
large  number  of  taps  are  needed  in  order  to  reject  the  input  noise. 

Eventually  algorithm  noise  becomes  the  significant  factor. 

2.  The  figures  show  only  the  mean  values  of  the  steady-state  weights  at  a 
particular  instant  of  time  after  the  filter  has  converged.  Hence,  there 
should  be  a quadratic  phase  shift  with  the  tap  number  in  accordance  with 

D(n)  in  Eq.  (6).  Comparison  of  Figures  42-44  with  Figure  41  and  Figures  40-48 
with  Figure  45  shows  that  the  steady-state  weights  do  display  this  behavior. 

3.  The  filter  trades  off  coherent  Integration  (proportional  to  the  number  of 
significantly  non-zero  weights)  against  the  phase  changes  required  at  each 
Up  to  follow  the  chirped  signal.  Since  the  phase  change  required  at  each 
iteration  for  each  tap  grows  linearly  with  tap  number  (entries  in  V), 
weights  at  the  far  end  of  the  line  must  make  large  phase  changes  in  com- 
parison to  those  at  the  beginning  of  the  line.  Note  that  the  quadratic 
phase  correction  along  the  line,  D(o),  is  independent  of  time.  Hence, 
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once  the  filter  estimates  fQ  and  f,  it  knows  D(o)  and  can  introduce  these 
phase  corrections  statically.  On  the  other  hand,  the  filter  must  change 
phase  by  the  entries  in  V at  each  iteration.  Large  phase  changes  are 
most  easily  made  when  the  magnitude  of  the  weights  are  small.  In  Fig- 
ures 42-44,  <10  = u>52  = 2 /10M  radians  and  in  Figures  46-48,  <10=  2 /40  M 
radians.  Hence  the  weight  has  to  change  by  it/ 5 and  t/20  radians, 
respectively.  In  order  to  accommodate  these  large  phase  changes  for  the 
same  algorithm  step  size,  the  weights  of  the  far  end  of  the  line  must  be 
smaller  than  those  at  the  beginning  of  the  line.  As  f decreases,  the  dif- 
ference in  phase  changes  at  the  two  ends  of  the  line  decreases  and  the  filter 
can  make  use  of  significant  values  for  the  weights  at  the  far  end  of  the  line. 
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THE  MEAN  SQUARE  ERROR  IN  STEADY-STATE 


T 

The  error,  € (n),  is  the  difference  between  d(n)  and  the  filter  output,  W (n)  X(m.  Its 
mean  square  value  is  given  by 


E|[d(n)  - WT(n)  X(n)]  [d*(n)  - W+(n)  X*(n>]  j 

E[d(n)  d(n)*]-  EjwT(n)  X(n)  d*  (n)]  - E[d(n)  W+(n)  X*(n)] 

+ E[wT(n)  X(n)  W+(n)  X*(n)] 


Using  the  assumptions  preceding  Eq.  (2), 


E[|e(n)|2]  = (cr.2  + a/)-  o2  2 Re  *E  [W(n)]+  D(n)  j + E[wT(n)  X(n)  W+(n>  X*(ni]  (28) 


The  middle  term  in  brackets  in  Eq.  (28)  can  be  evaluated  using  Eq.  (25)  and  is  given 
by 


E[WssH+  D(n) 


_ 2 Y'/  j JL  k62  ...  2.\ 

*°s  fev  ‘ n ' 

Ti  : - 


(29) 


The  last  term  in  Eq.  (28)  can  be  evaluated  as  follows.  Let  the  weight  vector  be 
written  as  a mean  value  plus  a zero- mean  fluctuation  process. 


W(n)  *E(W(n)]  + ;(n) 


(30) 


E[wT(n)  X(n)  W+(n)  X*(n)]  = E[w+(n)j  e[x*(d)  XT(n)]  e[w(ii)] 

+ E[;+(n)  X*(n)  XT  (n)  ;(n)] 

The  first  term  in  Eq.  (31)  is  known.  The  second  term  in  Eq.  (31)  is 
E[r+(n)  X*(n)  XT(n)  f(n)]  * e[?  +(n)  E [x*(n)  XT(n)]  c (n)] 

85  (°-s2+an2)E[f+(n,J(n)] 


* (as2+(7n2)Maw2  * 


assuming  the  weight  fluctuations  are  stationary,  uncorrelated  from  tap-to-tap,  and 

2 

have  the  same  variance  <rw  for  each  individual  weight.  Thus,  using  Eqs.  (31)  and 
(32)  in  Eq.  (28)  yields 


E jt(n)|2  = og2  1 - E|W(n)|  + D(n)  2 + ctr2  [l  + E|w(n)]+E|w(n) 


+ «,J2  . _ 2 


M (cr  2 + o 2 
\ s n , 


2 - (M+l)  pa 


w here  a has  been  approximated  by  the  weight  fluctuations  under  noise  alone 5’  ' ^ 
The  first  term  in  Eq.  (33)  represents  the  error  in  estimating  the  chirped  complex 
exponential  signal.  The  second  term  is  the  sum  of  the  noise  power  in  the  reference 
channel  and  the  noise  power  passed  by  the  mean  weights  of  the  adaptive  filter.  The 
last  term  represents  the  weight  misadjustment  variance  multiplied  by  the  total  input 
power. 

<y 

Eq.  (33)  normalized  by  the  total  input  power,  has  been  evaluated  as  a function  of  paQ" 
for  M = 16,  32,  64,  signal-to-noise  ratios  of  0 and  +10  dB  and  various  t2.  In 
Figures  49-54,  the  trade-off  can  be  seen  between  static  and  dynamic  contributions  to 
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total  mean-square  error.  In  each  case,  there  is  an  optimum  selection  of 
which  minimizes  the  deleterious  effects  of  signal  errors,  input  noise  and  weight 
misadjustment  noise.  Comparison  of  the  filter  performance  for  increasing  input 
signal-to-noise  ratios  verifies  improved  system  performance.  On  the  other  hand,  for 
sufficiently  large  dif  2 and  SNR  = 10  dB,  it  is  seen  that  the  normalized  man  square 
grror  increases  as  M increases.  This  effect  is  due  to  the  weight  mis— adjustment 
noise  exceeding  the  longer  coherent  integration  gain  obtained  with  larger  filters.  As 
it2  decreases,  a point  is  reached  where  sufficient  smoothing  time  is  available 
(small  (i  a 2)  to  reduce  the  weight  misadjustment  noise  to  a level  so  that  improved 
performance  is  obtained  for  longer  filters  (e.g.  w ? = 4 x 10"5). 

9 

It  can  be  seen  from  Figures  49-54  that  the  optimum  selection  of  nFn  , for  a given  filter 
length  M and  signal-to-noise  ratio,  varies  in  the  same  manner  asu  ( . As  !•  t 
increases,  a larger  value  ofpan2  is  required  to  achieve  the  minimum  mean-square 
errorr  However  this  minimum  mean-square  error  increases  as  a.  . The  filter 
has  less  time  to  learn  the  statistics  of  the  signal  and  hence  must  make  a larger 
mean-square-error  as  the  price  for  tracking  a faster  moving  signal.  It  is  somewhat 
difficult  to  determine  the  point  where  the  LMS  algorithm  loses  track  of  the  signal 
since  mean-square  error  is  not  a good  measure  of  tracking  performance. 
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CONCLUSIONS 


A mathematical  model  of  the  mean  weight  behavior  for  the  LMS  adaptive  filter  has 
been  presented  when  the  filter  is  operating  as  a single  frequency  line  enhancer  and 
line  follower.  For  a fixed  frequency  complex  sine  wave  input,  the  LMS  filter  weights 
have  been  shown  to  respond  to  signal  and  noise  more  rapidly  than  to  noise  alone.  This 
implies  that  the  filter  learns  more  quickly  that  a line  has  appeared  than  it  is  able  to 
forget  that  the  line  is  turned  off. 

When  the  signal  frequency  is  changing  linearly  with  time,  the  mathematical  model 
predicts  a time-varying  behavior  of  the  filter  mean  weights  necessary  to  respond  to 
the  changing  signal  frequency.  As  the  chirp  rate  increases,  the  filter  reduces  the 
relative  amplitudes  of  the  weights  so  as  to  adjust  the  effective  filter  length  to  optimally 
match  the  properties  of  the  signal.  That  is,  for  example,  suppose  the  filter  designer 
selects  a filter  of  length  M = Mj.  However,  the  chirp  rate  is  sufficiently  large  so  that 
the  chagge  ir.  signal  frequency,  between  algorithm  iterations,  is  greater  than  the  band- 
width of  the  filter.  Then,  the  LMS  algorithm  will  automatically  scale  the  amplitude  of 
its  weights  to  have  an  effective  length  Mg,  Mg  <M^,  such  that  the  signal  remains 
inside  the  adaptive  filter  between  iterations.  As  long  as  the  signal  frequency  lies 
within  the  LMS  filter  bandwidth,  the  LMS  filter  algorithm  can  track  the  changing  fre- 
quency since  there  is  sufficient  correlation  between  the  two  inputs  to  drive  the  LMS 
algorithm  in  the  correct  direction. 

The  mathematical  model  of  the  mean  weight  behavior  has  been  used  for  selecting  p, 
the  adaptation  coefficient  of  the  algorithm,  for  a wide  variety  of  signal  and  noise 
parameters.  The  criteria  of  optimality  was  that  of  minimizing  the  filter  output  mean 
square  error,  since  the  error  is  the  driving  term  in  the  weight  adjustment  algorithm. 
(An  alternate  criteria,  based  on  a signal  detection  model  using  the  filter  output,  could 
also  be  a candidate  for  optimization.)  A set  of  curves  of  normalized  mean-square 
error  as  a function  of  signal- to-noise  ratio  and  chirp  rate  were  obtained.  From  these 
curves,  the  following  observations  can  be  made: 

9 

1.  For  a given  signal-to-noise  ratio  and  chirp  rate  u><5“,  there  exists  an 
optimum  selection  of  p that  minimizes  the  overall  mean  square  error. 

2.  For  slowly  changing  signal  frequency,  the  nvan  square  error  exhibits  a 
relatively  broad  minimum.  This  is  because  a large  range  of  pctq  will 
follow  the  slowly  changing  signal  frequency  yet  allow  sufficient  smoothing 


u 
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On  the  other  hand,  for  a rapidly  moving  signal  frequency,  mismatch  in 

A 

selection  of  ^cr  can  cause  a significant  increase  in  mean-square  error  as 

B 2 
compared  to  the  optimum  selection  of  |j.crQ  . 

3.  As  a system  designer,  one  would  choose  a ^crQ2  that  would  be  optimum  for 


the  fastest  chirp  rate  expected.  The  mean  square  error  would  always  be 


upperbounded  by  the  mean  square  error  for  the  fastest  chirp  rate. 
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APPENDIX.  EIGENVALUES  OF  TRANSFORMATION  MATRIX 


In  order  to  easily  evaluate  the  01th  power  of  a matrix,  the  eigenvalues  of  the  matrix 
are  needed.  The  matrix  in  brackets  in  Eq.  (15)  is 

Po  v.  p^1  [1  - ^x]  . Po  V.  [p0->  - P/1 »] 

* P0  V*  ^1  - ^(Tn2  M(0)j  P0-1  (I-l> 

Since  PQ  premultiplies  and  PQ_1  post-multiplles  R=V*  [i-^2  M(0)j,  it  is  only 
sary  to  find  the  eigenvalues,  A,  of  R.  The  eigenvalues  of  R satisfy: 


neces- 


[i  - X*  k2  1 °oO] 


A I 


(1  - ,on2)  I - AV  - per/  DoDo+ 


- 0 


(1-2) 


where 


Dq  * D(0)  and  ] denotes  the  determinant. 


Because  of  the  simple  structure  to  R,  an  expression  for  the  eigenvalues  can  be  found. 
Given  a matrix  of  the  form  B * A + where  a and  b are  column  vectors, 


B|  = |A I [l  + bj*  A-1  aj 


(1-3) 
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Thus,  with  A * (1  - n.an2)  I - AV,  ^ • DQ  ■ b1# 


1 


j1  - »as2 


m 


a 


M 

Ji  *m 


-a) 


(1-4) 


The  term  in  brackets  yields  an  M'th  order  polynomial  in  A for  which  there  is  no 

general  analytic  solution.  Eq.  (1-4)  must  be  programmed  on  a digital  computer  for 
2 2 

various  p.,  ag  , on  and  a. 

Although  explicit  values  of  A are  not  obtainable,  a simple  upper  bound  on  the  eigen- 
values of  A and  hence  on  the  transient  behavior  of  Eq.  (21)  are  obtainable.  This 
upper  bound  on  the  eigenvalues  is  useful  since  it  is  an  Indication  of  the  slowest  possible 

response  of  the  system. 

Let  u be  an  eigenvector  of  the  matrix  in  Eq.  (1-1)  with  associated  eigenvalue.  A j. 
Then,  with  u + u = 1,  and  uT  = fujU^  . . . umj 

■Ss*  F.T'V1  (‘-’“'n2')  a tM) 

Now, 

* + . _ 
a i h aAi  - 

(I-pon2X)  u (I’6' 


U (I-p<Tn  \) 


-1 


p v*  V 

o o 


Using  PQ  ^ * PQ+  and  V * V 


2,  |2 


N2  - 1 !%l2l1->“’/xi 

j=l 


"l|2  1 - *(’ n * M<,s2)  + l“'n2|  S Kl 

j=2 


Let  q 


= ^ I uj  2 S 1 and  p = 1-Q  = 


<1.  Then 


2 I 2 

At|2  = <l-q)jl-  H (^n2+M<Ts2)|  + |l-  ^<Tn2| 

■ |»-^("-„2*M<Ts2)i’  i’q(1'l“,n2|2  - |1->1(<’n2*M,’82)i  ) 
s ! l- i“’n2i t<>r 


Similarly 


A[2  ■ P I1’  (<,n2,'M<’»2)(2  * ll"p)  I1'  l“n2|2  S I*"  “ (<’n2*“"*a)| 

for|l-n(i>n2»M<rs2)|>  l-pun2 

Eqs.  (1-8)  and  (1-9)  lead  to  the  following  bounds: 

(i)  If  p (on2  + M<ts2)  £ 1,  assuming  crn  t 0 

|Ai|2s|1-^(Tn2|<1 


(I-10> 
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(§)  If  k(<7d2  ♦ Mag2)>lbut  n<7n2  £ 1,  then  for  <rQ2  / 0 
| A J 2 * | 1-  p.<rn2| 2 < 1 for  S 2(1-^) 


and 


2 

< It*  (*n2+M<782)|  for  HM<7s2  > 2 (1-Han2) 

In  the  latter  case,  |aJ2<  1 if  n^<7n2+Mcrs2)  < 2 
(3)  If ^^n2>  1 *en 

li|2  S|1"‘i(an2+Mv)i 


(1-11) 


(1-12) 


and 


!<  1 if  ^(u^+MOg2)  < 2 

Combining  these,  It  can  be  seen  that  A.I  < 1 if  ^ (a  2 + Mcr  2 ) < 2.  Since  I A , J < 1 is 

M \ n s ' / 2 1 2\ 

the  condition  for  existence  of  a steady  state  mean  weight  vector,  p.  ^crn  +Mcs  f < 2 Is 
a sufficient  condition  for  a steady  state  solution.  It  is  interesting  to  note  that  the  bound 
in  each  case  is  just  the  magnitude  of  the  largest  eigenvalue  In  the  stationary 


on 


‘il 


2 2 

frequency  case,  with  a , a , n and  M unchanged.  Further 

Q S 


<2 


is  the  condition  for  the  convergence  of  the  adaptive  canceller  In  the  fixed  sinusoid 
case.  This  leads  to  the  surprising  conclusion  that  if  the  mean  weight  vector  achieves 
a steady  state  value  in  the  stationary  case  for  given  crg  , oQ  , M and  p.,  then  a steady 
state  solution  will  exist  for  those  parameter  values  regardless  of  the  rate  of  change 
of  frequency. 
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A more  direct  result  for  convergence  can  be  obtained  as  shown  below.  This  latter 
approach  does  not  bound  the  eigenvalues  however. 


From  the  discussion  following  Eq.  (1-1)  It  is  seen  that  R is  similar  to  Q and  therefore 

fl81 

has  the  same  eigenvalues.  From  1 \ 


lm  y Am_1  = U-  A f 1 

-•0O  ^ 


11m 

m=l 


(1-13) 


if  the  L_  norm  of  A is  less  than  unity,  l.e. , A <1.  In  our  case,  R = A and 


iRl  - V* 


[w^2  (X)] 


(1-14) 


But  II V* II  = 1 and  ||l-u.cr  2 a||  can  be  evaluated  explicitly  as  the  square  root  of  the 

^ n r ' 2 V T 2 1 

largest  eigenvalue  of  the  matrix  I - (iffn  A 1 1 - (icn  Aj. 

But  jl  - aJ  is  self-adjoint  and  has  only  two  distinct  eigenvalues. 

Ai = 1 +Mos2) 


A2  = A3  " XM  = 


(1-15) 


Hence 


1 “ ^fa2  + Mos“  ^ <J> 


and 


1 " 


<1 
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which  implies  j A.  ^ J < 1 if 

0<^(<rn2  +m<t82)<2 

in  agreement  with  the  discussion  following  Eq.  (1-12). 
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appendix  vrn 

COMPARISON  OF  ADAPTIVE  TRACKER  PERFORMANCE 
WITH  AN  EXISTING  AUTOMATIC  TRACKER 

It  is  of  interest  to  compare  the  performance  of  the  adaptive  tracker  with  a non-adaptive 
tracking  system.  For  this  purpose,  the  automatic  tracking  system  described  in 
reference  PJ  will  be  used.  Figure  55  shows  a block  diagram  of  the  tracker  to  be 
considered.  The  input  signals  from  the  left  and  right  half  beams  are  hard  clipped 
and  sampled,  with  theclipped  signal  for  one  side  subjected  to  avariable  delay,  r.  This 
delayed  signal  and  the  other  clipped  input  are  processed  by  a two  point  correlator, 
as  shown,  which  computes  the  correlation,  0 (t),  of  the  clipped  signal  at  r-  T.  The 
value  of  9 (T)  is  either  +1  or  -1  for  the  clipped  inputs.  The  contents  of  the  delay 
register,  which  determines  the  value  of  the  adjustable  delay,  is  increased  by  At  when 
+1  occurs  on  input  2 and  -1  on  input  1.  It  is  decreased  by  At  when  +1  occurs  on  input 
1 and  -1  on  input  2,  with  no  change  in  count  when  the  two  agree.  The  tracker  is  in 
steady  state  when  9 ("-T)  and  0 (t+T)  are  equal  (both  +1  or  both  -1).  Assuming  that  the 
correlation  function  of  the  input  is  symmetrical,  this  means  that  the  value  of  t is  equal 
to  the  delay  between  the  two  split  beam  inputs. 

The  performance  of  this  tracker  has  been  analyzed  in  reference  ClJ . The  analyzer 
is  in  excellent  agreement  with  both  computer  simulations  and  at-sea  tests  of  the 
device.  For  a broadband  stationary  target,  the  variance  of  the  delay  estimate  was 
shown  to  be 


where 


AT  = adjustable  delay  step  size 

cr  2 = signal  power 
s 

2 

cr  » noise  power 
n 

» 

0 __(T)  = value  of  the  derivative  of  the  normalized  autocorrelation  function 

99 

of  the  signal  at  r = T 


VIII- 1 


_ 


and  T = fixed  bulk  delay  used  In  the  correlator. 


This  resuit  is  valid  when  the  value  of  a ^ is  small  in  comparison  to  the  width  of 
*ss(T)* 

The  structure  of  this  tracker  is  considerably  different  than  that  of  the  adaptive 
tracker.  Typically,  the  delay  step  size,  at,  is  set  much  finer  than  the  time  between 
taps  in  the  adaptive  filter.  This  very  high  sample  rate  is  tolerable,  however,  because 
the  system  uses  one  bit  representations  of  the  input,  so  hardware  is  minimal.  The 
adaptive  tracker  replaces  a high  sampling  rate  with  multibit  processing  and  inter- 
polation between  coarser  delays.  In  spite  of  the  differences,  the  comparison  of  the 
two  is  interesting  in  light  of  the  wide  application  of  trackers  similar  to  the  clipped 
system. 


In  order  to  compare  realizations  of  the  two  tracking  schemes  that  are  in  some  sense 
equivalent,  the  time  constants  of  the  two  will  be  made  equal.  For  low  signal-to-noise 
ratio  the  time  constant  for  convergence  of  the  weight  vector  of  the  adaptive  filter, 
both  in  mean  and  variance  is  in  seconds 

V - * O]'1  Tat  «' 

where  T t is  the  adaptive  tracker  sample  interval. 

Reference  Cl  3 has  shown  the  time  constant  of  the  clipped  tracker  of  Figure  55  to  be 


t =. 
ct 


tr 

S 


(ir/2) 
5 — 


s n 


where  Fct  is  the  clipped  tracker  sampling  rate. 


(3) 
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It  Is  now  necessary  to  set  or  such  that  r = t^.  In  practice,  the  input  to  the  adaptive 
tracker  is  normalized  such  that 


This  normalization  will  also  be  applied  to  the  input  of  the  clipped  tracker.  Then, 
setting  (2)  equal  to  (3)  yields 


A r =. 


“5OTT- 


FctTat  <Tfi§R><-*'ss<T>> 


and  substituting  in  (1)  gives 


SNR  2 2 

16  pctVrria> 


where  SNR  = <r  2/<r  2. 

s n 


(4) 


(5> 


This  is  the  variance  of  the  time  delay  estimate  of  a clipped  tracker  whose  time 
constant  has  been  set  equal  to  the  adaptive  tracker,  and  is  not  valid  at  high  SNR. 


In  part  VI  of  this  report,  the  performance  of  an  adaptive  tracker  for  a broadband 
signal  with  an  800  Hz  flat  low  pass  spectrum  was  considered.  Following  the  design 
procedure  for  the  clipped  tracker  given  in  reference  Cl  1 gives  the  parameters 

Fct  = y = 2400  Hz 
|®ss*T)|  = 728.15  sec/sec 

For  the  adaptive  tracker  considered,  T t = (1/2400).  Substituting  these  parameters 
in  (5)  gives 

= (1. 1634  x 10“6)  M t^j|P)2  (G) 
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It  should  be  stressed  that  the  expressions  (5)  and  (6)  are  not  a generalized  compari- 
son of  the  particular  clipped  tracker  with  the  adaptive  tracker.  The  comparison  is 
between  an  adaptive  filter  and  a clipped  tracker  whose  time  constant  has  been 
adjusted  at  each  signal-to-noise  ratio  to  be  equal  to  that  of  the  adaptive  tracker. 

This  means,  roughly,  that  the  two  trackers  should  be  able  to  track  dynamic  targets 
with  the  same  rate  of  change,  or  that  the  time  to  achieve  the  predicted  bearing 
accuracy  for  an  emerging  target  should  be  the  same.  Given  the  great  difference  in 
the  structure  of  the  two  trackers,  however,  this  equivalence  is  only  approximate,  and 
should  be  used  only  for  order  of  magnitude  comparisons,  and  not  to  predict  fine 
differences  in  sensitivity. 

Reference  p.3:  C.  N.  Pryor,  "A  Simplified  Automatic  Tracking  Technique  for 

Signal  Correlation  Systems,"  NOL  Technical  Report  No.  67-152, 

21  September  1967. 
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